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Abstract 



In this thesis we investigate quantum mechanical effects to various aspects of gravita- 
tional collapse. These quantum mechanical effects are implemented in the context of the 
Functional Schrodinger formalism. The Functional Schrodinger formalism allows us to 
investigate the time-dependent evolutions of the quantum mechanical effects, which is 
beyond the scope of the usual methods used to investigate the quantum mechanical cor- 
rections of gravitational collapse. Utilizing the time-dependent nature of the Functional 
Schrodinger formalism, we study the quantization of a spherically symmetric domain 
wall from the view point of an asymptotic and infalling observer, in the absence of ra- 
diation. To build a more realistic picture, we then study the time-dependent nature of 
the induced radiation during the collapse using a semi-classical approach. Using the do- 
main wall and the induced radiation, we then study the time-dependent evolution of the 
entropy of the domain wall. Finally we make some remarks about the possible inclusion 
of backreaction into the system. 
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1. Introduction 



This thesis is based on work done in a series of nine papers, which have been pubhshed in 
several different journals, see Refs. [Ml Ell ESI EH 113 US El ES] • The work here is varied 
and takes several different aspects into account, however, for this thesis we concentrate 
on only a subset of these papers. The subset of interest here are those papers which 
include gravitational collapse of a massive shell only. Even though this is only a subset 
of the possible parameters that a black hole can have, this subset displays most of the 
interesting features that illustrate the core of our work. 
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2. Formalism 



The most important part of our research is the formahsm used to study the quantum 
mechanical effects of gravitational collapse. To study these quantum effects of gravita- 
tional collapse we will institute the Functional Schrodinger equation. In this part of the 
thesis we will first derive the Functional Schrodinger equation, which will be the primary 
equation used to study the quantum effects of gravitational collapse. 

The main purpose of the Functional Schodinger equation is to introduce the "observer" 
time into the Wheeler-de Witt equation. It is well known in General Relativity that for 
different foliations of space-time, different physical observations occur. For example, one 
can consider the equations of motion for a black hole. Let us consider an object which 
is falling into a black hole from different points of view. If one chooses the time as 
observed by an asymptotic observer as the desired foliation of space-time, upon solving 
the equations of motion of the object, one finds that it takes an infinite amount of time 
for the object to fall into the black (even if the object is a photon). However, if one 
choses the time as observed by a freely falling observer (one that is falling into the black 
hole along a geodesic) as the desired foliation of space-time, upon solving the equations 
of motion of the object in this case, one finds that it takes a finite amount of time for 
the object to fall into the black hole. Therefore, one can see that it is important and 
instructive to consider different foliations of the space-time to learn different aspects of 
the system of gravitational collapse. The Functional Schrodinger equation allows for one 
to specify the particular foliation of space-time that is of interest and study the system 
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CHAPTER 2. FORMALISM 



from that view point. 

A second purpose of the Functional Schrodinger equation is to allow one to investigate 
the time evolution of the system. Typically this is not done in the study of gravitational 
collapse. The preferred method of study is to consider an initial static asymptotically flat 
space-time, let the system evolve (with no knowledge of the evolution), then consider 
a final different static asymptotically flat space-time. Then by comparing these two 
different space-times, one can in principle have some understanding of what the evolution 
was like between these two events. The Functional Schrodinger equation will in principle 
allow us to study the total time evolution of the system, not just the static asymptotically 
flat regions of space-time. 

2.1. Functional Schrodinger Equation 

In this section we will derive the Functional Schrodinger equation. 

The Wheeler-de Witt equation for a closed universe is given by, see Ref.[T], 



where H is the total Hamiltonian and ^ is the total wavefunction for all the ingredients 
of the system, including the observer's degrees of freedom denoted by O. Eq. (j2.ip is a 
consequence of the idea that there is no "God" time, or no preferred time, or no super- 
observer time. Therefore Eq. ()2.ip is written in a gauge independent fashion, since there 
is no preferred observer to observe the system. 

In general we can write the wavefunction in Eq. ()2.1|) as 



Here X°' = X"{(^^) describes the location of the wall as a function of the internal wall 



H-^ = 



(2.1) 



M' = ^(X",5^„$,0). 



(2.2) 
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world volume coordinates (""i 9iiu is the metric, and $ is a scalar field. The Roman 
indices go over the internal domain wall world volume coordinates and the Greek indices 
go over space-time coordinates. Note that the wavefunctional in Eq. ()2.2p is a functional 
of the fields but not the space-time coordinates. In general, the total Hamiltonian is a 
linear combination of the Hamiltonian of the system itself and that of the Hamiltonian 
of the observer. Therefore we will separate the Hamiltonian into two parts, one for the 
system and the other for the observer, which can be written as 

H = Hsys + Hobs- (2-3) 

Any (weak) interaction terms between the observer and the wall-metric-scalar system 
are included in Hsys- The observer is assumed to not significantly affect the evolution of 
the system and vice versa. In mathematical language this means that we are assuming 
that the Hamiltonian for the system and the observer commute with each other 

[Hsys, Hobs] =0- (2.4) 

The total wavefunction Eq. ()2.2p can be written as a sum over eigenstates 

^ = Y,Ck^sysisyS,mLiO,t) (2.5) 
k 

where k labels the eigenstates and are complex coefficients. 

To solve the full Wheeler-de Witt equation is very difficult since it involves all the 
degrees of freedom, both that of the system and the observer. Here we shall utilize 
the frequently employed strategy of truncating the field degrees of freedom to a finite 
subset, hence we will be consider with the minsuperspace version of the Wheeler-de Witt 
equation. As long as we keep all the relevant degrees of freedom that are of interest, 
this is a useful truncation. Since we are only considering a subset of the total degrees of 
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freedom, this is now considered an "open" system. In an "open" system, one can then 
define an appropriate "observer" time in which one chooses to make measurements. 
Therefore we can write the Schrodinger equation for the observation as 

^ (2.6) 

This is convenient, however, we wish to make observations on the system not on the 
observer. To transform this to observation on the system we wih make use Eq. (|2.ip . 

To introduce the observer time on observations of the system, we use Eq. (j2.ip . Eq. (j2.5p 
and Eq. (j2.6p . therefore we can write 

k 

= Hsys CA:^',s(s2/^, t) + Hobs Ck^%s{syS. O^Ll^, 

k k 

k k 



where we made use of Eq. (|2.3p and in the last Une we used Eq. (j2.6p . 
Now consider the integral of the last term in Eq. (|2.7p , we have 



k 



dt Y Ck'fsysisyS, m^LiO, t)=Y Ck'fsysisyS, t)^obsiO, t) 



k 



k 

(2.8) 

However, by virtue of the Wheeler-de Witt equation, the total wavefunction ^ is time- 
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independent. Therefore the first term on the right hand side in Eq. (j2.8p is zero. Shrinking 
the integral we then have, 



k k 

Substituting Eq. (j2.9p into Eq. (j2.7p we can then write, 

Hsys^Ck^%Asys,t)^',,,{0,t) = iY,Ck (dt^%si^ys,t)) f (2.10) 



or interms of one k value, we then arrive at the Functional Schrodinger equation 

Hsys^sys = (2-11) 

For convenience, from now on we will denote the system wavefunction simply by ^ and 
drop the superscript k and the subscript "sys". Similarly H will now denote Hgyg, and 
the Schrodinger equation reads 

ff* = i^. (2.12) 

2.2. Discussion 

Here we have derived the Functional Schrodinger equation. As discussed above, the 
purpose of the Functional Schrodinger equation is to introduce the "observer" time into 
the Wheeler-de Witt equation, Eq. (|2.ip . This will allow us to be able to use the classical 
Hamiltonian of the system of gravitational collapse, then study the evolution of the 
system from the view point of any observer of our choosing. This has two benefits: First 
the formalism allows us to choose the "observer" we wish to study. As discussed earlier, 
different "observers" will observer different phenomena which are of interest. Secondly, 
the formalism will allow us to evolve the system quantum mechanically over time. One 
of the benefits of this approach is that we can, in principle, observe thermodynamic 
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properties of the system in a time-dependent fashion, which we will discuss in Chapters 
[7] and [HI As we will discuss, this is something which is beyond the scope of the usual 
methods used to study the thermodynamic properties of the system. 
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3. Occupation Number 



Throughout the text we will be interested in the number of particles created during 
the gravitational collapse of our object, i.e. the radiation. Therefore we will derive the 
occupation number of the particles created during the time of collapse in this chapter 
for future convenience. Throughout this text we will be interested in systems with 
spherical symmetry, since this is the simplest case to consider. Here we note that due 
to the spherical symmetry, gravitational radiation is excluded from the system, thus 
the radiation which we will consider will be from the excitation of particles due to the 
time-dependent nature of the gravitational metric. 

To consider the radiation we will consider a quantum scalar field (j) in the background 
of the gravitational collapsing object, which is given by 



where S is the action of the scalar field. The reason we are considering only a scalar 
field is that this is the simplest and easiest case, which gives insight into most of the 
physically significant phenomena. By considering more complicated fields, one arrives 
at the so-called gray-body factors, see for example Ref.[23], which are dependent on 
the type of field used. Here we derive the number of particles induced as a function of 
observer time "t" . Here "t" is used for any foliation of space-time used throughout this 
body of work, whether the time is that of an asymptotic observer or that of an infalling 
observer. 




(3.1) 
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In most cases we arrive at the Hamiltonian of the system from Eq. ()3.ip . which is of 
the form of a sum of uncoupled simple harmonic oscillator. To simplify the notation, we 
consider one eigenmode of the simple harmonic oscillator given by 



2m 2 ^ ^ 



(3.2) 



where p is the momentum conjugate to x and x is the eigenmode. Using the standard 
quantization procedure, upon inserting Eq. ()3.2p into Eq. ()2.12p . we can then write 



1 9^ 2/ N 2 



ip{x,t) 



' dt ■ 



(3.3) 



Here Eq. ()3.3p can be solved exactly by utilizing the invariant operator method first 
developed by Lewis and Reisenfeld, see Ref.[5] and Appendix |X1 Using this method, 
Dantas, Pedrosa and Baseia showed, see Ref.[3], that the exact solution to Eq. ()3.3p at 
late times is given by 



1/4 



exp 



(3.4) 



where pt denotes the derivative of p{t) with respect to t, and p is given by the real 
solution of the non-linear auxilarly equation 



1 



ptt + UJ {t)p = 



(3.5) 



with intitial conditions 



Pt{ti) = 



where ti is the initial time. The time-dependent phase a is given by 



(3.6) 



a{t) 



1 /■* dt' 



2Ju PHt'Y 



(3.7) 
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To find the occupation number of the induced radiation, consider an observer with 
detectors that are designed to register particles of different frequencies for the free scalar 
field <I> at earlier times. Such an observer will interpret the wavefunction of a given mode 
X at late times in terms of simple harmonic oscillator states, {fn}, at final frequency a). 
Here Co, is the value of the frequency evaluated at a time tf as seen by the observer. The 
number of quanta in eigenmode x can be evaluated by decomposing the wavefunction 
Eq. (j3.4p in terms of the states {^Pn}, and by evaluating the occupation number of that 
mode. To implement this, we start by writing the wavefunction for a given mode at time 
t > tf in terms of the simple harmonic oscillator basis at t = 

^p{x, i) = Cn{t)ipnix) (3.8) 
n 

where 

Cn = y dxLp*^{x)il}{x,t) (3.9) 

is the overlap, i.e. inner product, between the initial and final state of the wavefunc- 
tion. The occupation number at eigenfrequency oj by the time t > tf, is given by the 
expectation value 

N{t,oj) = ^n\cn\^ . (3.10) 

n 

To evaluate the sum in Eq. ()3.10|) . we use the simple harmonic oscillator basis states 
but at a frequency a) to keep track of the different u^s in the calculation. To evaluate the 
occupation numbers at time t > tf, we need only set Co = to{tf). So the simple harmonic 
oscillator basis states are written as (see for example Appendix A. 4 of Ref. |22j) 




(3.11) 
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where Tin are the Hermite polynomials. Then Eq. (j3.9p and Eq. ()3.4p together gives 



1 



1/4 gta 



1 



"n! J 



1/4 



where 



(3.12) 



(3.13) 



To find /„ consider the corresponding integral over the generating function for the 
Hermite polynomials 



J{z) 



die 




-22(l_2/P) 



Since 



we can then write 



£ — ^ n\ 



n=0 



2 = 



Therefore 



Since 



27r / 2\"/2 , , 




^„(0) = (-l)"/2V^l!l__Il!:, ,j = e^,en 

vn! 

and ^n(O) = for n = odd, we find the coefficients c„ for even values of n, 

(_l)n/2g^a |y / 2\"/2 (n- 1)!! 

Cn = — A/ — I 1 



(cDp2)l/4 V P 



P 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



(3.19) 
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For odd n, Cn = 0. 

We can now find the number of particles produced during the cohapse. Let 



X 



2 



(3.20) 



Then using Eq. (|3.10p we have 



N{t,u>) 



E 

n=even 

2 



n \Cn\ 



2 



d \ - 

X— 2^ 

n=even 
d 1 

X 



(n- 1)!! 



nv. 



-X' 



^/Of^\P\ dx ^\ - x^ 



X 



V^|P| (l-x2)3/2' 



(3.21) 



Now inserting Ea. (fXT3|) and Eq.dS^Q]) leads to 



1 \ 2 / \ 2' 



(3.22) 



3.1. Discussion 



Here we derived the occupation number of the radiation induced during the time of grav- 
itational collapse. The occupation number is measured by an observer with a detector 
at late times t > tf. As stated earlier, this was done for convenience since we will use 
this quantity several times during this text. 
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4. Model 



To study a concrete realization of black hole formation we consider a spherically sym- 
metric Nambu-Goto domain wall (representing a shell of matter) that is collapsing. To 
include the possibility of (spherically symmetric) radiation, as discussed in the previous 
chapter (Chapter [3|), we consider a massless scalar field, <I>, that is coupled to the grav- 
itational field but not directly to the domain wall. The action for the system is then 
given as 



The first term is the Einstein-Hilbert action for the gravitational field, the second is the 
scalar field action, the third is the domain wall action in terms of the wall world volume 
coordinates, (a = 0, 1, 2), the wall tension a, and the induced world volume metric 



As stated in Chapter [21 the coordinates = X^i^C,"") describe the location of the wall. 
The term Sobs in Eq. (|4.ip denotes the action for the observer. 

As discussed earlier, a general treatment of full Wheeler-de Witt equation, Eq. (j2.ip . 
is very difficult. So, we shall use the frequently employed strategy of truncating the 
field degrees of freedom to a finite set, typically including only the relevant degrees of 
freedom. In other words, we will consider the minisuperspace version of the Wheeler-de 




(4.1) 



lah = g^^.daXf'dtX 



(4.2) 
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Witt equation. As long as we keep all the relevant degrees of freedom, this is a useful 
truncation. Since we are considering spherically symmetric domain walls, we will assume 
spherical symmetry for all the fields. Thus, the wall is described by the radial degree of 
freedom R{t) only. 

The metric for the wall is then taken to be the solution to Einstein equations for a 
spherical domain wall. In Ref.[l] the metric, as follows from the spherical symmetry, 
outside the wall is given by 

ds^ = - (^1 - dt^ + (^1 - 'dr2 + rW, r>R{t) (4.3) 

where Rg = 2GM is the Schwarzschild radius in terms of the mass M of the wall, and 

dn"^ = dO"^ + sm'^9d(l)'^. (4.4) 

By Birkhoff 's theorem, in the interior of the spherical domain wall the line element of 
the metric is flat, i.e. Minkowski, which is given by 

ds"^ = -dT"^ + dr^ + r^dO^, r < RH). (4.5) 

Here T is the interior time coordinate, not to be confused with temperature. The interior 
time coordinate is related to the asymptotic observer time coordinate t via the proper 
time r of the domain wall. By matching the coordinates for the interior and exterior at 
the wall, in analogy with the Isreal junction condition (see Ref.[5]), and assuming that 
the wall is infinitely thin, we have the relations 



dT / /di?^ ^ 
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and 



where 



dt 1 /„ /di?^^ 



B s 1 - |l. (4.8) 

By taking the ratio of Eq. (|7.42p and Eq. ()4.7p . the relationship between the interior time 
T and the asymptotic time t is given by 



dt -/bTM V B 

where Rr = dR/dr and R = dR/dt. 

Since we are restricting the system to fields with spherical symmetry only, we need 
not include other metric degrees of freedom. Thus, the scalar field can also be truncated 
to be the spherically symmetric modes 

$ = $(r,t). (4.10) 

In Ref.[l], Ipser and Sikivie integrated the equations of motion for the spherically 
symmetric domain wall. They found that the mass is actually a constant of motion and 
is given by, see Appendix|B]for a sketch of the method used, 

1 



M = -AttgR^ 
2 



yJl + Rl + yjB + i?2 (4.11) 



where it is assumed that max(i?)> I/AttctG. This assumption is just used to ensure that 
one does not start off inside of the collapsing spherical domain wall. 

By virtue of Eq. (j4.8p . Eq. (|4.1ip is implicit since Rs = 2GM. Solving for M explicitly 
in terms of Rr gives 



M = A-KaR^ 



Vl + -R? - 27raGR . (4.12) 
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However, making use of Eq. (j4.9p we can solve for M in terms of Rt = dR/dT 



M = AnaR^ 



1 



2TTaGR 



(4.13) 



1 — R^ 



Before we proceed we wish to discuss the physical relevance of Eq. ()4.12p . First consider 
the case where Rr = 0, i.e. for a static domain wall. The first term in the square bracket 
is just the total rest mass of the shell. When the shell is moving, i.e. Rr ^ 0, the first 
term in the square bracket takes the kinetic energy of the domain wall into account. The 
last term in the square bracket is the self-gravity, or the binding energy of the domain 
wall. Therefore we can identify Eq. (j4.12p (Eq. ()4.13p ) as the total energy of the system, 
hence the Hamiltonian of the system. Thus, we will refer to Eq. ()4.12p (Eq. ()4.13p ) as the 
Hamiltonian. 

4.1. Discussion 

Here we developed the classical Hamiltonian for a massive spherically symmetric domain 
wall undergoing gravitational collapse. As stated in the last paragraph, Eq. ()4.12p is the 
conserved mass of the system, however, it can be interpreted as the Hamiltonian. For 
the remainder of the text, we will use Eq. (|4.12p (Eq. ()4.13p ) as the Hamiltonian for the 
system. 
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In this chapter we wish to study the classical equations of motion of the collapsing 
spherically symmetric domain wall. To do so, we will consider the cases for two different 
foliations of space-time. Here we consider the collapse of the spherically symmetric 
domain wall from the point of view of an asymptotic observer (one who is at rest with 
respect to the collapse) and from the point of view of an infalling observer (one who is 
riding along with the shell). This will be done by considering Eq. ()4.13p and Eq. ()4.12p . 

A naive approach to obtaining the dynamics for the spherical domain wall is to insert 
Eq. (|4.3p and Eq. (|4.5p . as well as Eq. (|4.12p into the original action Eq. ()4.ip . Upon doing 
so it is known that this approach does not give the correct dynamics for gravitating 
systems. Therefore we will take an alternative approach for finding the action. We will 
find the action that does in fact lead to the correct mass conservation law. The form of 
the action can be deduced from Eq. ()4.13p fEq. ()4.12p ). 

5.1. Asymptotic Observer 

First we will consider the equations of motion from the view point of the asymptotic 
observer. The asymptotic observer is any observer stationary with respect to the col- 
lapsing domain wall, typically taken to be located at infinity. Here we summarize the 
work originally done in Ref . |30j . 

From Eq. (|4.13p we find the effective action for the spherically symmetric domain wall 
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to be 



l-R^- 2TTaGR 



(5.1) 



Using Eq. ()4.9p we can write Eq. (|5.ip in terms of the asymptotic observer time t as 



S, 



iira j dtR^ 



(5.2) 



From Eq. (j5.2|) the effective Lagrangian for the system is 



A-KaR^ 



(5.3) 



The generahzed momentum 11^ can be derived from Eq . (15.31) in the usual manner, 
this is given by 



AnaR^R 



1 



27rcrGi?(l - B) 
Vb^-r"^ Jb'^-{1-B)R? 



(5.4) 



Therefore from the Lagrangian, Eq. ()5.3p . and the generahzed momentum, Eq. ()5.4p . the 
Hamiltonian can then be written as 



H = iiraB^/'^R^ 



1 



2TTaGR{l - B) 
- i?2 Jb^-{1-B)R^ 



(5.5) 



For later convenience, we wish to find the Hamiltonian as a function of (R, Hr). To do 
so, we need to eliminate R in favor of Hr using Eq. (|5.4p . This can be done, in principle, 
but is very messy (the solution will involve solutions of a quartic polynomial). However, 
we will be interested in what is happening as the shell approaches the horizon, i.e. when 
R is close to Rs and hence B ^ 0, since this is the most interesting region of study. In 
this limit one can see that the denominators of the two terms in Eq. (|5.4p are equal. So 
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we can rewrite Eq. (|5.4p as 



where 

fi = a{l- 2TTaGRs). (5.7) 



In the region ~ i?^ the Hamiltonian, Eq. ()5.5p . is then approximately given by 

H . (5.8) 

= V(^n^)' + ^(4^/^^')'- (5.9) 

Here we note that the Hamiltonian, written in the form of Eq. ()5.9p . has the form of the 
energy of a relativistic particle with a position dependent mass. 

Since the mass is a constant of motion, the Hamiltonian is a conserved quantity, so 
from Eq. (|5.8p we can write 

B^^R^ , , 

h = , (5.10) 

\/s2 - m 

where h = H/Att^ is a constant. 
Solving Eq. ()5.10p for R we obtain 



R = ±B\ll-^^. (5.11) 



In the region R ^ Rg, this takes the form 



Since in the region i? ~ i? — )■ 0, the dynamics for the collapsing spherically 
symmetric domain wall in this region can be obtained by solving the expression 

R = ±B. (5.13) 
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To leading order in R — Rg, the solution is 

R{t) ^Rs + {Ro - Rs)e^''^' (5.14) 

where Rq is the radius of the shell at t = 0. Since we are interested in the collapsing 
shell, we take the negative sign in the exponential term, Eq. (|5.14p . 

Here we wish to make some comments on Eg . (15.141) . By virtue of the negative sign in 
the exponential, this then implies that, from the classical point of view, the asymptotic 
observer never sees the formation of the horizon of the black hole, since Eq. (j5.14p equals 
Rg only as t — 7- oo. This is in agreement with the fact that it takes an infinite amount 
of time for a photon to reach the horizon of a pre-existing black hole, as seen by an 
asymptotic observer (see for example Ref.[7]). Therefore, Eq. (|5.14p makes sense from 
this point of view. 

In Figure [5T] we plot the position of the domain wall for the asymptotic observer. Here 
we see the asymptotic behavior of the time dependence of the position of the domain 
wall. As shown in Eq. (j5.14p . the domain wall asymptotes to the Schwarzschild radius, 
taking an infinite amount of time for the domain wall to reach the Schwarzschild radius. 




^ i/Rs 

2 4 6 8 



Figure 5.1.: Here we plot the the solution in Eq. ()5.14p . 
Figure 15.21 shows the corresponding velocity of the domain wall as seen by the asymp- 
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totic observer. Here we see that the velocity of the domain wall asymptotes to zero as 
the domain wall collapses toward the Schwarzschild radius, as given in Eq. ()5.13p . 

m/dt)/Rs 




Figure 5.2.: Here we plot the the solution in Eq. ()5.13p . 



5.2. Infalling Observer 

Now we turn our attention to the infalling observer case, where the conserved mass is 
given by Eq. (j4.12p . Here we point out the misnomer in the name infalling. The infalling 
observe here is not to be confused with the traditional view point of an infalling observer, 
one who is traveling along a geodesic, or a freely falling observer. The observer in this 
case is infalling from the fact that the observer is attached to the domain wall and is 
infalling with the wall. Therefore, eventhough the observer is in a locally Minkowski 
reference frame, the overall reference frame is still Schwarzschild, since at any point in 
time the observer is in a Schwarzschild reference frame. Here we summarize the work 
originally done in Ref. [31 j. 

The effective action consistent with Eq. ()4.12p is 

Seff = -^-rra J drR^ y/TTR^ - RrSmh-'^{Rr) - 2-KaGR . (5.15) 

Therefore the effective Lagrangian expressed in terms of the infalling observer's time r 
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is given by 



Leff = -AiraR^ ^/l + i?^ _ smli'^Rr) - 2-KaGR 



From Eq. ()5.16p the generalized momentum 11^ is derived to be 



(5.16) 



lip, = AiiaR^ sinh"^(i?^) 



(5.17) 



From Eq. (|5.16p and Eq. ()5.17p . the Hamiltonian in terms of Rr is given by 



H = AnaR^ 



y/l + R^ - 27raGR 



(5.18) 



which is just Eq. (j4.12p as expected. 
From Eq. (|5.18p we can calculate R^ 



Rr = ± 



\ 



h 



^2 + 2^'^^^ j 



(5.19) 



where h = H/Ana. In general, Eq. (|5.19p cannot be solved analytically, at least not in 
very nice way. However, we can take some special cases to investigate the behavior of 
the solution to find the time dependence. 

As a first case we consider the zeroth order behavior near the horizon, i.e. R ^ Rg. 
In this region we note that we can write Eq. (|5.19p as 



Rr = ± 



\ 



+ 2TraGRs 



(5.20) 



hence Rr is constant. Integrating Eq. (|5.20p gives 



R{t) = Ro-t^ 



+ 2iTaGRs 



(5.21) 
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where Rq is the radius of the sheh at r = 0. 

As a second case we consider the case that h/R'^ » 2-KaGR > > 1 . Therefore we can 
write Eq. ()5.19p as 



R - ^ 



Integrating Eq. ()5.22p we then have the solution 



(5.22) 



R{t) = [R^- 3hT 



1/3 



(5.23) 



Eq. ()5.23p then gives that the time for an infaUing observer to reach Rg is 



Rl-Rl 



3h 



(5.24) 



Here we make some comments on Eq. (|5.2ip and Eq. (j5.23p . These solutions imply 
that the infalling observer will reach Rs in a finite amount of his/her proper time. This 
result is expected from classical general relativity, since the observer is in a locally flat 
Minkowski reference frame. Therefore, there is no difficulty for the observer once he/she 
reaches the horizon, the horizon is just another locally flat point in space according to 
this observer. 




Figure 5.3.: Here we plot the numerical solution to Eq. (|5.19p . 
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For consistency, in Figure 15.31 we plot the numerical solution of Eq. (j5.19p for the 
parameters h = 1/2 and a = O.lRs- Figure [531 shows that the observer does in fact 
reach Rs in a finite amount of his proper time. Figure 15.41 compares the special cases 
discussed above with that of the numerical solution. Here the blue curve is the full 
solution of Eq. ()5.19p . the green curve is the solution to Eq. ()5.20p and the red curve is 
the solution to Eq. ffO^]) . 

R/Rs 




12 3 4 

Figure 5.4.: Here we plot a comparison of the different approximations for the solution 
of Eq. (j5.19p . The blue curve is the solution to Eq. (|5.19p , the green curve is 
the solution to Eq. ()5.20p and the red curve is the solution to Eq. ()5.22p . 

In Figure 15.51 we plot the numerical solution for the velocity of the domain wall as 
seen by the infalling observer, from Eq. (|5.19p . Since the domain wall crosses its own 
Schwarzschild radius at a time of r = 1.66, Figure [531 shows that the velocity is infact 
approximately constant as R ^ Rs- After the domain wall crosses the Schwarzschild 
radius, the velocity then increases and diverges as i? — t- (the classical singularity). 

In Figure 15.61 we plot the numerical solution for the acceleration of the domain wall 
as seen by the infalling observer, from Eq. (j5.19p . Since the domain wall crosses its own 
Schwarzschild radius at a time of r = 1.66, Figure [5T6l shows that the acceleration is 
increasing almost linearly as R ^ Rs- After the domain wall crosses the Schwarzschild 
radius, the acceleration then increases and diverges as i? — )■ (the classical singularity). 
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dR/dr 




T 



Figure 5.5.: Here we plot the corresponding numerical solution for Eq. ()5.19p . Here we see 
that as R ^ Rs (t = 1.66), the velocity of the domain wall is approximately 
constant. However, after the domain wall passes the Schwarzschild radius 
the velocity diverges as i? — 0. 

Therefore, we can conclude that even though the observer is attached to the domain wall, 
he/she is not a truly free- falling observer (as stated at the beginning of this Section). 
The infalling Schwarzschild observer is an accelerated observer during the entire duration 
of the collapse. 



5.3. Comparing Asymptotic versus Infalling 

Here we wish to investigate the discrepancy of the observation between the infalling and 
asymptotic observers. 

To understand this discrepancy we turn to Eq. ()4.7p . As discussed in the previous 
section, in the limit R ^ Rg, Eq. (|5.19p is approximately constant. This means that we 
can then write 

At « ^VB + cAt. (5.25) 
Now, since i? — )• in this limit and the proper time taken to reach the Schwarzschild 
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Figure 5.6.: Here we plot the corresponding numerical solution for the acceleration as- 
sociated with Eq. (j5.19p . Here we see that as R ^ Rs {t = 1.66), the 
acceleration of the domain wall increases almost linearly. However, after 
the domain wall passes the Schwarzschild radius the acceleration diverges 
as i? ^ 0. 



radius is finite, we can then see that 



lim At oo. (5.26) 

R^Rs 

This is a fairly crude approximating, thus in Figure \57l\ we plot dt/dtau versus R/Rs- 
Figure [5771 shows that as R — )• Rg, dt/dr does indeed diverge as given in Eq. (j5.26|) . 
Therefore Eq. (|5.26p can be thought of as the gravitational red-shift, which is the source 
of the discrepancy between the observations between the two observers. In Figure ISTSl we 
show the the position of the domain wall as a function of time for both the asymptotic 
and infalling observers. Figure ESI shows that initially the two observers are in agreement 
on where the domain wall is compared to the Schwarzschild radius. However, as the time 
increases (both asymptotic, t, and infalling, r) the discrepancy becomes more apparent. 
According to the asymptotic observer, the domain wall asymptotes to the Schwarzschild 
radius, while according to the infalling observer this happens in a finite amount of time 
(as stated earlier). 
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dt/dT 




Figure 5.7.: Here we plot dt/dr versus R/Rg- Here we can see that as i? — > i?s, dt/dr ^ 
oo as given in Eq. (|5.26p . 




0.5 1.0 1.5 2.0 2.5 3.0 



Figure 5.8.: Here we plot a comparison of the position of the domain wall relative to 
the Schwarzschild radius for both of the two different observers. Here the 
position of the domain wall as seen by the infalling observer is given in blue, 
while that of the asymptotic observer is given in green. 
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5.4. Discussion 

In this section we investigated the classical equations of motion for the collapsing spher- 
ically symmetric domain wall. As discussed in the first section, the asymptotic observer 
sees the domain wall collapse to the Schwarzschild radius Rs only as t — )■ oo. This is 
not an unreasonable result since in classical General Relativity an asymptotic observer 
never sees a photon cross the Schwarzschild radius since it take an infinite amount of 
time t to reach the horizon. Therefore one can easily believe the result here, since the 
time taken is due to the gravitational redshift of the photon. As the shell approaches 
Rs the redshift will increase until it becomes infinite by the time it reaches the horizon. 

In the second section, the infalling observer will see the shell collapse to the horizon 
in a finite amount of time. As discussed earlier, this is because the infalling observer is 
always in a locally flat Minkowski frame. Thus, the horizon is not a significant point for 
the observer, therefore there is no problem for him/her to pass right through and not 
even know it. This result is again expected from classical General Relativity. 
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In this chapter we wish to study the quantum equations of motion using the Functional 
Schrodinger equation Eq. (|2.12|) . again for the two different fohations of space-time. The 
idea here is that the quantum mechanical effects will change some of the difficulties that 
arise from the classical solutions. The ultimate goal is that examining these quantum 
mechanical effects will give us insight into the quantum mechanical nature of gravita- 
tional collapse and will possibly help guide us to be able to construct the appropriate 
theory of quantum gravity. 

Some of the difficulties that arise from the classical solutions of gravitation collapse 
are discussed below. This is not meant to be an exhaustive list, but give the reader of 
an idea of the topics that we wish to address in this text. 

First we will discuss one difficulty faced by the asymptotic observer, the presence 
of the horizon. Why is the event horizon a difficult place for the asymptotic observer, 
while there is it is no problem for the infalling observer, since the horizon is nothing more 
than a coordinate singularity? Under the classical notion, as discussed in Chapter [5l the 
presence of the coordinate singular creates an apparent gravitational time dilation, which 
makes the collapsing domain wall appear to stop. This effect is due to the divergence of 
the coordinate singularity when the observer sees the collapsing domain wall reach the 
Schwarzschild radius. There has been much discussion in this matter about how these 
quantum mechanical corrections can eliminate this effect. One such idea that has gained 
much attention over the past few years, we will discuss this process here. The basic idea 
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here is that the upon quantizing the sheh, the shell will now have quantum fluctuations 
in the position of the horizon. These fluctuations will then imply that the position of the 
horizon is no longer fixed, but will now be given by Rg + 6Rs, where 6Rs represents the 
small fluctuations in the position of the horizon. These effects can then make the time 
as measured by the asymptotic observer finite (see for example Sec. 10.1.5 of Ref.[8]) 



If this were correct, we would be able to observe black hole formation (due to the collapse) 
and other effects in finite time. Note however, that the fluctuations can go either way. 
In the case of Rg — 5Rs, the result becomes infinite agin. 

Now we will discuss difficulties facing the infalling observer. For this observe there are 
really two important regimes: The regime in the region R Rs and the region ~ 0. 
For the region ii ~ i?^, the concern is the exact opposite of that of the asymptotic 
observer. Will the quantization of the shell contradict the classical observation that the 
infalling observer sees the shell collapse to Rg in a finite amount of time? For the region 
ii ~ 0, classically this represents the point of the classical singularity. Penrose and 
Hawking showed in Ref.[9] that singularities are endemic in classical General Relativity. 
The question then arises whether these singularities are an intrinsic property of space- 
time or simply reflect our lack of the ultimate non-singular theory. The general belief is 
that quantization will rid gravitation of singularities. This is analogous to another theory 
where quantization got rid of the singularity, Electromagnetism. In atomic physics the 
singularity of the Coulomb potential, which has an identical 1/r behavior, was eliminated 
via quantization (see for example \10\ [TH [T2t [T3] ) . 

In this chapter we investigate these ideas. In later chapters we will investigate quantum 
mechanical corrects to addition aspects of gravitational collapse, those being thermody- 
namic. 




(6.1) 
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6.1. Asymptotic Observer 

First we will consider quantization of the shell from the view point of the asymptotic 
observer. This will be done by using Eq. ()2.12p . For this section we outline the work 
done in Ref . [30] . 

To utilize Eq. ()2.12p we need to use the Hamiltonian in Eq. ()5.9p . However, we notice 
that Eq. ()5.9p has a squareroot in it. Therefore, we will consider the Hamiltonian squared 

H'^ = BUrBUr + BiATTfiR^f. (6.2) 



Before we proceed, we discuss the choice of ordering in the first term on the right 
hand side of Eq. ()6.2p . Since we are considering quantum mechanics, the distance R 
and the conjugate momentum are now promoted to operators, which obey the standard 
commutation relations. Thus in general we would need to add terms to the squared 
Hamiltonian in Eq. ()6.2p that depend on the commutator [S,n/j]. However, in region of 
interest we find that the commutator is given by 

its 

Estimating H by the mass M of the domain wall, as for the discussion in Chapter HI the 
terms due to the operater order ambiguity will be negligible provided 

M » 

Rs M 

where mp is the Planck mass. Therefore we can ignore the ordering ambiguity and chose 
the ordering given in Eq. ()6.2p . provided that this limit is satisfied. 
Now we apply the standard quantization procedure. 



[R,IiR]=i. 
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We substitute 



n« = -,|^ (6.3) 



into the squared Schrodinger equation 



Inserting Eq. (|6.2p into Eq. ()8.2.1|) we then obtain 



To find the wavefunction for the collapsing domain wall we need to solve Eq. (|6.5p . To 
solve Eq. ()6.5p in terms of R can be a formidable exercise in mathematics. However, we 
can simplify the matter by defining the tortoise coordinate 



u = R + Rsln 



We can then see that Eq. (|6.6p then gives 



Rx 



(6.6) 



Bn„ = (6.7) 



where we used Eq. ()6.3p . Using Eq. (|6.7p we can then rewrite Eq. (j6.5p as 

+ B(47r/i/?2)^^' = 0. 



We can now identify Eq. (j6.8p as just the massive wave equation in a Minkowski back- 
ground with a mass term that depends on the position of the domain wall. We need to 
now solve Eq. (j6.8p for u. To do so, we must first write the mass term in terms of n, 
rather than its present state of R. However, some care is needed since at i? = i?s, n is 
divergent, so we must take the appropriate branch. From Eq. ()6.6p we have that for the 
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region R G oo) maps onto u G (—00, 00) and R G (0, Rs) maps onto u G (0, —00). 

We are interested in the situation of a collapsing shell, hence the region R G {Rg, 00). 
Thus we are interested in u in the region u G (—00,00). We can solve Eq. ()6.8p for the 
entire region, however, we are mostly concerned with the effect when ~ i?<j. In the 
region R ^ Rg, the logarithm in Eq. ()6.6p dominates, so we can write 

~ i?, + i^,e"/■^^ (6.9) 



We look for wave-packet solutions propagating toward Rs, or in terms of u, u — t- —00. 
Thus from Eq. (j6.9|) we have 

B ~ e"/^= ^ 0. (6.10) 

This means that the last term in Eq. ()6.8p . the mass term, can be ignored in this region. 

In the region ~ i? s , the dynamics of the wave-packet is simply given by the free wave 
equation, where any function of light-cone coordinates (n it) is a solution. To make 
this explicit, we consider a Gaussian wave-packet propagating toward the Schwarzschild 
radius 

* = _^e-("+*)'/2^' (6.11) 
V2^s 



where s is some chosen width of the wave packet in the u coordinate. The width of the 
wave-packet remains fixed in the u-coordinate while it shrinks in the R coordinate via 
the relation dR = Bdu, as follows from Eq. (j6.6p . 

Let us consider some properties of Eq. (j6.1ip . First, we see that the wave-packet 
travels at the speed of light in the u coordinate. This is expected since the Schrodinger 
equation takes the form of a massless wave equation in Minkowski space. Further, in 
the u coordinate, the wave packet must travel out to u = —00 to get to the horizon, 
R = Rs- Thus we can conclude that the quantum domain wall does not collapse to Rs 
in a finite amount of asymptotic time. 

Therefore one can conclude that the quantum mechanical effect, i.e. the quantization 
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of the domain wall, does not smear out the presence of the coordinate singularity at 
the Schwarzschild radius. The asymptotic will not see the formation of the horizon in a 
finite amount of his/her time. Hence, the quantum solution does not alter the classical 
result found in Chapter [5l 



6.2. Infalling observer 

Now we consider quantization of the domain wall from the view point of the infalling 
observer. This will be done using Eq. (j4.12p as the Hamiltonian of the system. As 
discussed earlier, here we wish to solve the Functional Schrodinger equation in two 
different regions of interest. The first region is that near the Schwarzschild radius, Rg. 
The second is the region near the classical singularity, in the region i? ~ 0. For this 
section we summarize the work originally done in Ref. [31 j. 

6.2.1. Near Horizon 

The exact Hamiltonian in terms of Rr is again given by Eq. (j4.12p . From Eq. (|4.12p we 
again see the presence of a square-root. However, in this case it is not as easy to remedy 
this as it was in the case of the asymptotic observer, since upon squaring the Hamiltonian 
will not get rid of the square root. To simplify the analysis we will require that Rr is 
small. This is indeed a restriction to the special motion of the wall, since in general 
Rr can be large near Rs if the shell is falling from a very large distance. However, one 
may always choose initial conditions in such a way that the initial position of the shell 
R{t = 0) is very close to Rs- 

In the limit of small Rr, Eq. ()4.12p simplifies to 

H = 4TTaRl l + lRl-2TTaGRs . (6.12) 
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Then again in the same hmit, the conjugate momentum Eq. ()5.17p simphfies to 

Ur = AiraRlRr. (6.13) 

Ignoring the constant terms from the Hamihonian Eq. ()6.12p and using Eq. (|6.13p we can 
write the Hamiltonian as 

Using the standard quantization procedure, we substitute 

nH = -i|j (6.15) 

into the Schrodinger equation Eq. ()2.12p . which yields 

1 _ M 



(6.16) 



Investigating Eq. (j6.16p . we see that Eq. (j6.16p is just the Schrodinger equation for a 
freely propagating "particle" of mass A-kgR^, as one can expect from this approximation. 
Since Rg is only a finite distance away for an infalling observer we conclude that the 
wavefunction will collapse at Rs in a finite amount of proper time. 

For the question on if the quantization of the domain wall will cause problems for the 
infalling observer, we can conclude that quantum effects do not alter the classical result. 
Hence a collapsing shell crosses its own Schwarzschild radius in a finite proper time. 

6.2.2. Near the Origin 

Now we wish to investigate the quantization of the domain wall in the region of the 
classical singularity i? ~ 0. 

The exact Hamiltonian in terms of Rr is again given by Eq. (j4.12p . where R-r is given 
by Eq. ()5.19p . In the region near the classical singularity, i.e. in the limit i? — t- 0, the 
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classical expression for Rr (keeping only the leading order term) becomes 

Rr - (6.17) 

where h is defined in Chapter [5l Eq. ()6.17p clearly shows that in this region the classical 
expression for Rr diverges. Up to the leading term near the origin, Eq. ()6.17p implies 
that the Hamiltonian is 

H = A-KcrR^Rr. (6.18) 

Substituting the asymptotic behavior, Eq. (|6.17p . in the expression for the generalized 
momentum Eq. (j5.17p we have 

liR = 4TTcrR^ sinh-^(i?^). (6.19) 

From Eq. (|6.19p we see that 

lim Ur = 0. (6.20) 



This then gives 

R-^b 47rcri?2 
This implies that R^^ which is defined as 



li"n7^ = -~- (621) 



by virtue of Eq. ()5.17p . near the horizon becomes 
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Therefore substituting Eq. ()6.23p into Eq. (j6.18p we have 



2'iTaR? exp 



( 



i d 
AnaR"^ 'dR 



^{R,t) 



' dr 



(6.24) 



Let us consider some properties of Eq. ()6.24p . The differential operator in the exponent 
gives some unusual properties to the equation. First we note the presence of the R^^ 
term. This implies that if we expand the exponent we can not stop the series after a 
finite number of terms, but instead need to include all orders of the expansion, we will 
make this explicit below. Thus, we need to include an infinite number of derivatives of 
the wavefunction ^ into the differential equation. An infinite number of derivatives of a 
certain function uniquely specifies the whole function. Thus, the value of (the derivative 
of) the function on the right hand side of Eq. ()6.24p at one point depends on the values 
of the function at different points on the left hand side of the same equation. This is in 
strong contrast with ordinary local differential equations where the value of the function 
and certain finite number of its derivatives are related at the same point of space. This 
indicates that Eq. (|6.24p describes physics which is not strictly local. 

Here we will make the non-locality of Eq. (|6.24p explicit. To illustrate we will use 
the expression involving the sinh term, this is not necessary, however, it will make the 
explanation more clear. Note that we can rewrite sinh as 



By making a change of variable we can make this more explicit. If we introduce the new 
variable v = R^, Eq. (j6.24p becomes 



We can then see that the the differential operator in the exponents in Eq. (|6.25p are 



sinh (a;) 



2 




(6.25) 
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just a translation operator, which shifts the argument of the wavefunction by a non- 
infinitesimal amount of 3i/47rcj. Since the wavefunction is complex in general, a shift by 
a complex value is not a problem. Therefore Eq. ()6.25p can be written as 



irav 



2/3 



^ ( U + ) - ^ ( 



3i 



.d^{v,T) 



(6.26) 



Here we make an additional change in variable and define v' = v — 2>i/ATTa. Then we can 
rewrite Eq. ()6.26p as 



/ 3i \ 

vrcr ( v + 

Ana 



2/3 



[v' + -^-.r ) - ^' {v' ,T 



M{v,t) 
dr 



(6.27) 



To interpret this we rely on usual differential calculus. From calculus we have 



J{x + Ax) - fix) ^ fix) + ^(Ax)"/W - /(a 

n=0 



^(Ax)"/(") « AxJ^. (6.28) 



Here the last step assumes that Ax is small so we are justified at keeping only the first 
term in the expansion. Now, performing the same procedure as in Eq. (|6.28p we can 
write 



ij {v' + Av' ,t) - {v' ,t) « '0(t.',T) + ^(At;')"V'^"^ -^(^^',t) 

n=0 

oo 

= ^(At>')"V^"^ (6.29) 



n=0 



where 



Au' 



6i 
47r(j 



(6.30) 



However, here we cannot truncate the series after a finite number of derivatives since 
by virtue of Eq. (|6.30p . Av' is not a small shift, provided that a does not go to infinity. 
Therefore we must keep all orders of the derivative, since as in Eq. ()6.28p each derivative 
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has higher powers of Av'. 

An interesting thing to note here is that as o" — t- 0, the non-local effect becomes 
stronger and more predominant. A possible understanding of this is as follows. Outside 
of the domain wall, there exists a certain Hilbert space, while on the inside there exists 
a second Hilbert space. The transition from the first Hilbert space to the second is not 
necessarily a smooth transition. When the domain wall collapses to the singularity, the 
effect of each of these Hilbert spaces is now taken into account. When dealing with 
a massive domain wall, the warping of space-time is greater, however it exists over a 
larger distance. This makes the transition more smooth from point to point. However, 
for a small domain wall, the warping is less noticeable at larger distances. Therefore, 
the transition is more violent the closer one gets to the classical singularity in this case, 
since in both cases the warping diverges. 

What about the value of the wavefunction at the classical singularity? Eq. (|6.26p shows 
that the wavefunction near the origin ^'(i? — t- 0, r) is in fact related to the wavefunction 
at some distant point ^'(i? — )• {^^)^^'^ ,t). This also implies that the wavefunction 
describing the collapsing domain wall is non-singular at the origin. Indeed, in the limit 
R ^ 0, this equation becomes 

d^{R^O,T) 

-q:^ = (6-31) 

where we used the fact that the wavefunction at some finite R, i.e. ^{R — )• (^^)^/^,t), 
is finite. From Eq. (|6.31|) it then follows that ^{R — )■ 0) = const. This gives strong 
indication that quantization of the domain wall may indeed rid gravity of the classical 
singularity. 



39 



CHAPTER 6. QUANTUM TREATMENT 



6.3. Discussion 

In this chapter we quantized the collapsing domain wah and investigated the quantum 
corrections. In the first section we did this with respect to the asymptotic observe to 
see if these fluctuations changed the classical observation that the domain wall takes an 
infinite amount of time to reach Rg. Upon quantizing the domain wall, we found that 
in this view point the classical scenario was not changed by the quantum fluctuations. 

In the second section we investigated the quantization of the domain wall from the 
view point of the infalling observer. Here we did this for two different points of interest, 
near the horizon i? ~ i?^ and near the classical singularity i? ~ 0, respectively. In the 
region R ^ Rg, we found that upon quantizing the domain wall, the classical view point 
was again unchanged. The quantum fluctuations did not alter the fact that according to 
the infalling observer, the domain wall will collapse to Rs in a finite amount of proper 
time. In the region ~ 0, we found some interesting properties of the wavefunction. 
First, we found that the physics of the wavefunction in this region are strongly non- 
local, meaning that the value of the wavefunction at i? ~ depends on the value of the 
wavefunction some distance away from the classical singularity. This situation has been 
previously suggested in the context of the information loss paradox (see for example 
[14: \ \15 \ [16]). As we pointed out in the section, what is interesting is that this non-local 
behavior becomes increasingly manifest in the limit that o" — )• 0, i.e. that the mass of 
the domain wall becomes very small. This may be a consequence of the non-separability 
of the Hilbert space between the outside of the domain wall and the singularity. In 
the massive domain wall scenario, the transition from the outside to the inside Hilbert 
space is smoother than in the case of the light domain wall. Secondly, we found that 
the wavefunction is in fact finite at the classical singularity, which implies that quantum 
fiuctuations may rid gravity of the classical singularity. 
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In this chapter we wish to investigate one of the thermodynamic properties of gravita- 
tional coUapse. Two of the most important thermodynamic properties of a black hole 
are the temperature (discussed in this chapter) and the entropy (discussed in Chapter 
[8j). As mentioned in Chapter [2l the benefit of the Functional Schrodinger equation is 
that this will allow us to evolve the system over time. This is in contrast with the usual 
method used to evaluate the thermodynamic properties of a black hole. 

The most widely used method of determining thermodynamic properties of a black 
hole is the so-called Bogolyubov method. The method here is as follows. One considers 
an initial asymptotically flat space-time, usually Minkowski, at the beginning of the 
gravitational collapse. The system is then allowed to evolve to a final asymptotically flat 
space-time, Schwarzschild in the context of a shell of matter only, with no knowledge of 
what happens in between. Then by matching the coefficients between these two space- 
times, the mismatch of these two vacua gives the number of produced particles. As 
mentioned, what happens in between the initial vacua and the final vacua is beyond the 
scope of the Bogolyubov method. 

Since the Functional Schrodinger equation allows one to find the time dependent 
wavefunction for the system, one can, in principle, ask the question of what happens 
during the evolution of the collapse. In this chapter we will investigate the time evolution 
of the radiation, in the form of the occupation developed in Chapter [31 during the time of 
gravitational collapse. The occupation number will then allow us to fit the temperature 
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of the radiation, and compare with that of the pre-formed black hole. 

As discussed in Chapter U we can consider the radiation given off during the col- 
lapse of the domain wall by considering a massless scalar field ^ that is coupled to the 
gravitational field. The action of the scalar field can then be written as in Eq. (j4.ip 



We decompose the (spherically symmetric) scalar field into a complete set of real basis 
functions denoted by {fkir)} 



The exact form of the function /fc(r) will not be important for us. We will, however, 
be interested in the wavefunction for the mode coefficients {ak{t)}. Note, again here we 
use "t" to be anytime coordinate of interest to us, not necessarily the asymptotic time. 

From Eq. (j4.5p and Eq. (|4.3|) we can see that the action for both the different foliations 
of space-time will consist of two parts, one from Eq. ()4.5p and the second from Eq. ()4.3p . In 
both foliations, the asymptotic observer and the infalling observer, we will be interested 
in the region R ^ Rs, therefore we will explicitly write out Eq. ()7.ip then take the limit 
to find the dominating contributions. 

From this action we can then find the Hamiltonian of the system using the usual 
methods. Substituting the Hamiltonian into Eq. ()2.12p we can find the time-dependent 
wavefunction of the system. As discussed in Chapter [3l the time-dependent wavefunction 
allows us to find the occupation number N, which will allow us to fit the temperature 
of the radiation (which we will describe below). 

In 1975 Hawking showed that for a pre-existing static black hole, the black hole will 
radiate its mass away, see Ref.[T7j. The radiation that is given off has a finite tempera- 
ture, as viewed by an asymptotic observer, which is known as the Hawking temperature. 
Therefore, it will be instructive for us to compare our late time result with Hawking's 




(7.1) 




(7.2) 



k 
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original calculation. 

For the infalling observer, the calculation is also instructive to give us an idea of 
what the region is like for this observer. Some unanswered questions are: since the 
temperature as measured by the asymptotic observer is finite, what is the temperature 
at the horizon for the local observer? If the temperature is infinite at the horizon, as one 
would expect since the temperature for the asymptotic observer is finite, will the infalling 
observer burn up before he reaches the horizon? To answer this question, we will use 
two different foliations of space time, that of Schwarzschild and Eddington-Finkelstein, 
respectively. 

7.1. Asymptotic Observer 



Here we consider the radiation as measured by the asymptotic observer. For this section 
we summarize the work originally done in Ref . [30| . 

As stated above, the action for the scalar field can be written in two parts 



S — Sin + Sout 



(7.3) 



where 



27r 



J dt J drr^ 



Sout = 2-K j dt drr 

'm 



f 

1 - Rs/r 



+11-^ 



(7.4) 
(7.5) 



where T is given in Eq. (|4.9p . which with Eq. ()5.1ip . gives 



T = B\ll + {l-B)^. 



(7.6) 



As mentioned above, we are interested in the R ^ Rs behavior of the action. As R ^ Rs, 
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we see that T ~ i? — t- 0. Therefore the kinetic term in Eq. (j7.4p diverges as {R — Rs)~^ 
in this Umit, while the kinetic term in Eq. (|7.5p diverges logarithmically. Therefore the 
divergence of the kinetic term in Eq. ()7.4p dominates over that of the divergence of the 
kinetic term in Eq. (j7.5p . The gradient term in Eq. ()7.4p vanishes in this limit, while the 
gradient term in Eq. ()7.5p becomes finite. Thus the gradient term in Eq. ()7.5p is dominant 
over that of the gradient term in Eq. (|7.4p in this limit. Hence the action can be written 
as 



S ^ 2tt dt 



1 



Rs 



drr'^{dt^f + / drr"^ ( 1 
-D /o JRs 



R. 



(7.7) 



where we have changed the limits of integration to Rg since this is the region of interest. 
This approximation is valid provided the contribution from r G (Rs, R{t)) to the integrals 
remains subdominant, and also the time variation introduced by the true integration 
limit can be ignored. 

Now, using Eq. ()7.2p we write the action in Eq. ()7.7p as 



S 



dt 



-^dk{t)Mkk'dk'it) + ^afc(t)Nfcfc/afc'(t) 



where M and N are matrices that are independent of R{t) and are given by 



rRs 

Mkk'=A7T drr^Ur)fk'{r) 
Jo 

roo / 

Nfcfc/ =Att drr^ | 1 
JRs 



f'kir)fUr)- 



(7.9) 
(7.10) 



Using the standard quantization procedure and Eq. ()2.12p . the wavefunction il){ak-,t) 
satisfies 

.dip 



dt 



where 



. d 

I 

dak{t) 



(7.11) 



(7.12) 
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is the momentum operator conjugate to afc(t). 

The problem of radiation from the collapsing domain wall is equivalent to the prob- 
lem of an infinite set of uncoupled harmonic oscillators whose masses go to infinity with 
time. We can see from Eq. (j7.9p and Eq. (j7.10p that the matrices M and N are hermitian. 
Therefore, it is possible to do a principal axis transformation to simultaneously diago- 
nalize M and N (see Sec. 6-2 of Ref. [18]) for example). Then for a single eigenmode, 
the Schrodinger equation takes the form 



R 2m dh^ 2 



^(M) = .^ (7.13) 



where m and K denote eigenvalues of M and N, and h is the eigenmode. 
Dividing Eq. (j7.13p through by i?, we can write in the standard form 



^(M)=^^ (7.14) 



where 



and 



ri= fdt(l-?^] (7.15) 







R 



where we have chosen to set r]{t = 0) = 0. 

From Eq. ()5.14p we can see that the classical late time behavior of the shell is given 
by 1 — Rs/R ~ exp(— t/i2s). For early times, the behavior depends on how the spherical 
domain wall was created and we are free to choose a behavior for R{t) that is convenient 
for calculations and interpretation. The most convenient case to use is a static begin- 
ning. This can be obtained if we artificially take the collapse to stop at some time, tj. 
Eventually we can then take — >• oo, as given by Chapter [5l We will then chose R to 
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be 

1 te(-oo,o) 



e-*/«% tG(0,t/) (7.17) 
Q-tf/Rs^ te{tf,oo). 



With the choice of initial static space-time of the domain wall, the initial vacuum state 
for the modes is the simple harmonic oscillator ground state, 

V^(6,r? = 0) = (^)'^%— o^V2. (7.18) 

The exact solution for late times is given by Eq. (|3.4|) with initial conditions given by 
Eq.dMl). 

As discussed in Chapter [3] an observer with a detector will interpret the wavefunction 
of a given mode b at late times in terms of simple harmonic oscillator states at the final 
frequency 

id = uJoe^f/^^' (7.19) 

where we have made use of Eq. (j7.17p . 

The number of quanta in eigenmode b can be evaluated from Eq . (|3.21|) . By calculating 
it can be checked that remains constant for t < and also t > tf. Hence all the 
particle production occurs for < i < and is a consequence of the gravitational 
collapse. 

Now we can take the limit — )• oo. In this limit, p remains finite but — >• — oo as 
t > tf ^ oo, provided ojq / (see Appendix O for details). However, we are interested 
in the behavior of A^ for fixed frequency, a). From Eq. (|7.19|) in this limit implies ojq — ?• 0. 
From the discussion in Appendix [Cl we also know that p — )• oo as wq — )• 0. Hence we 
find 

Ar(t,w)~-^~— tytf^oo. (7.20) 
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Therefore the occupation number at any frequency diverges in the infinite time hmit 
when backreaction is not taken into account. In Figure \77i\ we have plotted the occupa- 
tion number N versus t /Kg for various values of u>Rs . Figure 17.11 confirms the late time 
behavior of the time dependence of the occupation number, Eq. ()7.20p . 

ojRs =1,5, 10, 20 




Figure 7.1.: N versus t/Rs for various fixed values of ooRg. The curves are lower for 
higher ojRg- 

We have also numerically evaluated the spectrum of mode occupation numbers at any 
finite time and show the results in Figure 17.21 for several different values of t /Rg ■ Figure 
17.21 shows that as the asymptotic observer's time increases, the occupation number of 
larger values of ciji?^ increases. This is consistent with Eq. (|7.16|) . since as i — t- oo, oj — t- oo. 

To find the temperature of the radiation, we compare the curve in Figure [7^2] with the 
occupation numbers for the Planck distribution 

^H^) = (7.21) 

where /3 is the inverse temperature. We can see that the spectrum of occupation numbers 
is non-thermal. As an example, there is no singularity in at a; = at finite time. 
However, as t — t- oo, the peak at u = does diverge and the distribution becomes more 
and more thermal for these times. There are also oscillations in N. 



47 



CHAPTER 7. RADIATION 



t/Rs = 6, 7, 8 




Figure 7.2.: N versus uRg for various fixed values of t/Rg. The occupation number at 
any frequency grows as t/Rg increases. 

We now wish to fit the temperature of the radiation. From Eq. (j7.2ip . we can find the 
inverse temperature to be 

,3 = !^<l±i/M=r-'. (7.22) 

(jjRs 

In Figure [7^ we plot ln(l + 1/A^) versus ORg for various values of t/Rg. Here we see 
that for smaller values of t/Rg the spectrum for /3 is non-thermal. For example, Figure 
17.31 shows a thermal-like distribution for only small values of ORg for t/Rg = 2, while 
the larger values or u)Rg are not yet thermally induced. If one fitted the slope of /3 for 
this particular time, the only relevant region is that between < CoRg < 200. We can 
also see that the fluctuations of /3 are large. However, as t/Rg increases more and more 
values of u>Rg are thermally induced, hence one can fit a larger region. The fluctuations 
for larger values of t/Rg become much smaller, until they become almost completely 
non-existent. Another feature which occurs as t/Rg increases is that the slope of /3 goes 
to zero, which would imply that the temperature of the radiation in fact goes to infinity, 
not to a finite number as predicted by Hawking. 

However, from Eq. ()7.14p we see that the time derivative of the wavefunction on the 
right-hand side is with respect to r], not with respect to t, and uj is the mode frequency 
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t/Rs = 7, 8, 9 

Ln(l + 1/N) 




Figure 7.3.: ln(l + versus ujRg for various values of t/Rg. 

with respect to 77 as well. Eq. (j7.15p tells us that the frequency in t is (1 — Rs/R) times 
the frequency in so at the final time tf, this implies 

Jt) ^ ^-tflRs^ (7^23) 

where the superscript (t) on cj refers to the fact that this frequency is with respect to 
time t. However, since we are interested in the temperature in time t, we must also 
rescale the temperature in the same manner as the frequency. So the temperature seen 
by the asymptotic observer is 

T = e-*//^=/3-^(t/). (7.24) 

Fitting a thermal spectrum to the collapsed spectrum of Figure 17. 2^ as shown in Figure 
17.41 we obtain 

19 

T « ^ = 2.4r^, (7.25) 

where Tn = l/AnRg is the Hawking temperature. Since there is ambiguity in fitting 
the non-thermal spectrum by a thermal distribution, we can only say that the constant 
temperature, T, and the Hawking temperature are of comparable magnitude. 
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t/Rs = 8 

Ln(l + 1/N) 




Figure 7.4.: ln(l + versus QRs for t/Rg = 8. The dashed hne shows ln(l + 1/Np) 
versus CoRs where Np is a Planck distribution. The slope gives /3~^ and the 
temperature is given in Eq. (|7.25p . 

7.2. Infalling Observer 

Here we consider the radiation as measured by the infalling observer. To do so we will 
consider two different foliations of space time. As seen in Chapter [5l the acceleration of 
the Schwarzschild observer becomes divergent as the observer crosses the horizon. There- 
fore it is important to switch to another observer whose acceleration is no longer diver- 
gent upon crossing the horizon. For this observer, we will work in Eddington-Finkelstein 
coordinates. In this section we summarize the work originally done in Ref. |32] . 

First we will consider the Schwarzschild observer, after which we will consider the 
Eddington-Finkelstein observer. 
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7.2.1. Schwarzschild Coordinates 



The action can again be written in two parts, see Eq. ()7.3p . where 



Sin =2vr I dr / drr 

/•oo 

Sout =2vr I dr drr"^ 
Jr{t) 



1 



B 



(7.26) 



2 1 - Rjr ^ B V r J"^ ' ' 



(7.27) 



The most interesting things happen when the sheh approaches the Schwarzschild radius. 
From Eq. (|5.19|) we see that Rr is constant in the hmit when R ^ Rg. Therefore 
the kinetic term for Sin is roughly constant. The kinetic term in Sout goes to zero as 
ii — )• so the Sin kinetic term is dominant. Similarly the potential term in Sin goes 
to a constant while the potential term in Sout becomes very large, so the potential term 
in Sout dominates. Therefore in the region R ^ Rg we can write the action as 



1 



drr , 



(7.: 



where we have changed the limits of integration from -R(r) to Rg since this is the region 
of interest. 

Using the expansion in modes, Eq. ()7.2p . we can write the action as 



S= dr 



1 1 \R I 

o / „., =«fc(T)Afcfc/dfc/(r) + -^ak{T)Ckk'ak'{T) 

2 ^1 + i?2 2B 



(7.29) 



where d = da/dr, and A and C are matrices that are independent of R{t) and are given 
by 



rRs 

Akk' = 47r / drr^fk{r)fk>{r) 
Jo 

Ckk' = ^TrJ^drr' (l - fk{r)fk,{r). 



(7.30) 
(7.31) 
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Prom the action Eq. (j7.28p we can find the Hamiltonian, and according to the stan- 
dard quantization procedure, the wave function Tp{ak,T) must satisfy the Functional 
Schrodinger equation. We can write the Schrodinger equation as 



dip 



\R'T 



Ofc(T)Cfcfc/afc/(T) 



(7.32) 



where 



a 



(7.33) 



dakir) 

is the momentum operator conjugate to ak{T). 

Again, the problem of radiation from the collapsing domain wall for the infalling 
observer is equivalent to the problem of an infinite set of uncoupled harmonic oscillators 
with time dependent mass and frequency. Following the principal axis transformation 
used in the section above, the single eigenmode Schrodinger equation take the form 



^(6,t) 



.avK^r) 

' dr 



(7.34) 



where m and K denote eigenvalues of A and C, and h is the eigenmode. 
Re- writing Eq. ()7.34p in the standard form we obtain 



drj 



(7.35) 



where 



and 



K \Rr\ 



m 



2 



jdr'^/l + R^ 



(7.36) 



(7.37) 



where we defined uoq = K/m. 

To proceed further, we will use the classical background of the collapsing domain wall 
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Eq. (j5.2ip . The initial vacuum state for the modes is the simple harmonic oscillator 
ground state, 

V'(6,r, = 0) = (^)'^%— o^V2. (7.38) 

The exact solution for late times is given by Eq. p.4p with initial conditions given by 
Eq.dMl). 

As discussed in Chapter [3] an observer with a detector will interpret the wavefunction 
of a given mode 5 at late times in terms of simple harmonic oscillator states at the final 
frequency O. 

The number of quanta in eigenmode b can be evaluated from Eq. ()3.2ip . By calculating 
Nj- it can be checked that N remains constant for r < and also t > Tf. Hence all 
the particle production occurs for < r < rj and is a consequence of the gravitational 
collapse. 

Now we can take the limit rj — )• Tc. In this limit, p remains finite but — t- — oo as 
T > Tf ^ Tc, provided ojq 7^ (see Appendix|D] for details). However, we are interested 
in the behavior of for fixed frequency, to. From the discussion in Appendix |Dl we 
also know that p — t- cxo as — 0. Therefore the occupation number at any frequency 
diverges in the infinite time limit when backreaction is not taken into account. 

In Figure 17.51 we plot the occupation number of produced particles as a function of 
time (for several fixed frequencies ORs). The amount of proper time needed for the shell 
to reach Rg can be obtained by integrating Eq. (j5.19p . For a = O.Oli?"'^ this critical 
proper time is Tc = 7/3Rs- Figure 1731 shows that, as the infalling observer approaches 
Rs, the occupation number increases and diverges exactly at Rs- The same conclusion 
as found by analyzing the occupation number A^ as a function of p and pr (see Appendix 
Id]) . This is in agreement with what one would expect in the absence of backreaction. 
Hawking showed, see Ref.|17j. that the flux of particles at late times diverges for a fixed 
background, i.e. fixed mass of the object. Here, from Eq. (j4.12p . we are treating the mass 
of the domain wall as a constant of motion. This means that we keep adding energy 
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to the domain wall during the time of collapse, despite the the loss of mass due to the 
radiation. For the asymptotic observer it takes an infinite amount of his time for the 
domain wall to collapse to Kg, see Chapter [5] for discussion of this. However, this infinite 
time interval corresponds to a finite amount of time for the infalling observer's time. 
Thus, one may conclude that the infalling observer has to encounter the infinite number 
of particles produced during this finite amount of time before he reaches Kg. 



cdRs = 2, 4, 10 




Figure 7.5.: The occupation number as a function of proper time r/Rg for various 
fixed values of particle frequencies ooRg. The curves are lower for higher 
values of uiRg. The occupation number diverges as the infalling observer 
approaches Rg, which happens as t ^ Tc- 

We have also numerically evaluated the spectrum of mode occupation numbers at any 
finite time and show the results in Figure 17.61 for several different values of r/Rg. Figure 
17.61 shows that as the infalling observer time increases, the occupation number of larger 
values of coRg increases. 

To find the temperature of the radiation, we again compare the curve in Figure 17.61 
with the occupation numbers for the Planck distribution, which is given by Eq. (|7.2ip . 
where /3 is again the inverse temperature. We can see that the spectrum of occupation 
numbers is non-thermal. As an example, there is no singularity in at c^j = at finite 
time. However, as r — )• Tc, the peak at w = does diverge and the distribution becomes 
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t/Rs = 2.33, 2.333, 2.3333 




ojRs 

Figure 7.6.: The occupation number as a function of frequency CoRg for various fixed 
values of proper time r/Rg. The occupation number increases for larger 
values of t/Rs as t ^ Tc. 

more and more thermal for these times. There are also oscillations in N, which are not 
present in the Planck distribution. 

We now wish to fit the temperature of the radiation, however, from Eq. (|7.35p we see 
that the time derivative of the wavefunction on the right-hand side is with respect to r], 
not with respect to r, and uj is the mode frequency with respect to ij as well. Eq. (|7.37p 



tells us that the frequency in r is y^l + i?^ times the frequency in i]. However, recall 
from Chapter [5] Eq. (j5.19p tells us that as R ^ Rs, Rt is in fact a constant. Therefore, 
rj and r, for the case of the infalling observer, only differ by a constant amount. Hence, 
without loss of generality, we can ignore this shift by a constant amount, since the general 
features of the temperature will be the same. From Eq. (j7.2ip . we can find the inverse 
temperature to be 

^^H^±M£l^T-\ (7.39) 

LoRs 

In Figure \77l\ we fit a thermal spectrum to the collapsed spectrum of Figure [7l6l Several 
important features of the Hawking-like radiation can be taken from from this plot. First, 
the non-thermal features of the radiation are apparent. However, the departure from 
thermality (the fluctuations) are larger for earlier times, hence larger frequencies. This 
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observation was first argued in Ref. |30j . Second, as r —t- Tc and the infalling observer 
approaches Rg, the radiation becomes more and more thermal even at large frequencies. 
Third, at r = Tc, i.e. R = Rg, the radiation becomes purely thermal. At this point, 
the black hole is formed and the radiation becomes thermal, as known from various 
studies of quantum radiation from a pre-existing black hole. Finally, it is apparent that 
the slope of ln(l + versus uRg is decreasing as the infalling observer approaches 

Rg. Exactly at Rg, the slope of the curve is zero, indicating that the temperature of 
the radiation is infinite. This is not surprising since, as it is well known, the asymptotic 
observer in the nearly flat asymptotic region will register Hawking radiation with a finite 
temperature (see previous section). When the temperature is blue-shifted back to Rg, 
it clearly diverges. 

t/Rs = 2.333, 2.3333, 2.33333 

Ln(l + 1/N) 




Figure 7.7.: Plot of ln(l + as a function of frequency ORg for various fixed values 

of proper time t /Rg. The slope of the best fit line is /3, which is the inverse 
temperature. The non-thermal features disappear and the temperature di- 
verges as the Schwarzschild radius is approached, i.e. r — Tc- 

7.2.2. Infalling Eddington-Finkelstein Coordinates 

Now we consider the collapse from the point of view of an infalling Eddington-Finkelstein 
observer. This is a different space-time foliation than that in Schwarzschild coordinates, 
and we expect crucially different results. In particular, since the metric is not divergent 
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at the horizon, we do not expect infinite temperature there. 
For this purpose, we define the ingoing null coordinate v as 



v = t + r* (7.40) 



where r* is the tortoise coordinate. We can then rewrite Eq. ()4.3p as 



ds"^ = - {l- ^) dv^ + 2dvdr + r^dQ^, r > R{v). (7.41) 



where the trajectory of the collapsing wall is r = R{v). The interior metric is the same 
as in Eq. ()4.5p . The interior time coordinate, T, is related to the ingoing null coordinate, 
V, via the proper time on the shell, r. The relations are 



dT / /dRV 



and 



where 



dv I I dR I [dRV , ,^ 



B^l-^. (7.44) 



Consider again a massless scalar field $ which propagates in the background of the 
collapsing shell. The action for the scalar field is 

S = j d^x^^^g^'^d^'^d,^, (7.45) 

where g^^'^ is the background metric given by Eqs. ()4.5p and (|7.4ip . Decomposing the 
(spherically symmetric) scalar field into a complete set of real basis functions denoted 
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by {/fc(r)} 

^ = Y,ak{v)fk{r) (7.46) 

k 

we can find a complete set of independent eigenmodes {b^} for which the Hamiltonian 
is a sum of terms. 

Since the metric inside and outside of the sheh have different forms, we again spht the 
action into two parts 

Sin = 27T JdT Jj^'^ drr^ [-{dr^f + [dr^f] , (7.47) 



Sout = 271 I dv I drr^ 
lR(v) 



(7.48) 



We are again interested in the near horizon behavior of the radiation, i.e. as R ^ Rs- 
In this hmit we can write Eq. ()7.43p as 



dv 



1 



dr 2Rr 



(7.49) 



where Rr = dR/dr. Then with the help of Eq. (|7.42p we can write Eq. (|7.47p as 



Sin = 2tt I dv 



R(v) 



drr^ 



={dv^) 

2 ^RJ2{RJ2 + 1) 

+ 2^Rj2{R^/2 + l){dr^f 



(7.50) 



where Ry = dR/dv. Obviously, the action is not singular as R{v) — )■ unlike the 



58 



CHAPTER 7. RADIATION 



Schwarzschild case. From Eqs. (|7.48p and (|7.50p we can write the total action as 



S ^2tt / dv 



.1 



1 



drr'^ , 

2 y/Rj2{R„/2 + 1) 

oo roD 

+ I drr^d^<i>dr<i>+ / drr'^dr<^d.,<i> 

fOO / p 



+ / drr"^ ( 1 
JRs 



(7.51) 



where we have changed the hmits of integration from R{v) to Rg since this is the region 
of interest. 

Now using the expansion in modes Eq. ()7.46p . we can rewrite the action as 



dv 



1 



1 



2^R„/2{RJ2 + 1 



zOif^A-ff-i^idi^i 



(7.52) 



where d = da/dv, and A, Y and C are matrices that are independent of R{v) and are 
given by 



rRs 

Akk' = 27r drr'^fk{r)fk'{r), 
Jo 

/•oo 

Yfcfc, =47r/ drr^fk{r)fUr), 
JRs 

Rx 



Ckk' = Stt / drr^ 1 



Rs 



r 



f'k{r)f'k'ir)- 



(7.53) 
(7.54) 
(7.55) 



However if we take that the matrices are symmetric and real, we can see that Y = Y ^ , 
so we can write the action as 



S 



dv 



1 



1 



2 ^RJ2{RJ2 + 1) 



dk-^kk'dk' 



+ -^^kk' {dk^k' + «fc«A;') + -^ClkCkk'ak' 



(7.56) 
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Prom the action Eq. ()7.56p we can find the Hamiltonian, and according to the standard 
quantization procedure, the wave function •0(afc,f) must satisfy 



«^ = Hip, 
ov 



(7.57) 



or 



.dip 
' dv 



+ \k (Vi?,/2(i?,/2 + l)Y2^,(A-i)fcfe, + Ckk) at' 



+ -^/Rj2{Rj2TT)nkY 



ip 



(7.58) 



where 



Hi 



Oak 



(7.59) 



is the momentum operator conjugate to a^. Using the momentum operator conjugate 
to fl/fc, we can rewrite the Schrodinger equation as 



.dip 
dv 



^y/Rj2{Rj2 + l)nfc(A-i)fcfc/nfc, 



+ -ak { y^Rj2{Rj2 + l)Ylk'{A ^)kk' + Ckk' ) Qfc' 



i-V^,/2(i?,/2 + l)Yfcfc, ( )fcfc,(^fcfc, 



ip 



(7.60) 



where (5^^/ is the Kronecker delta function. 

So the problem of radiation from the collapsing domain wall for the infalling observer 
is equivalent to the problem of solving an infinite set of decoupled damped harmonic 
oscillators with time-dependent frequency. Since A, Y and C are symmetric and real, it 
is possible to simultaneously diagonalize them using the principle axis transformation. 
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Then for a single eigenmode, the Schrodinger equation takes the form 



.dip 
dv 



2 \ m J 



Zm 



(7.61) 



where m, y and K denote eigenvalues of A, Y and C, and b is the eigenmode. 
Re- writing Eq. ()7.6ip in the standard form we obtain 



Id m n . y 



ipib,r]) 



.dip{b, rf) 
drj 



(7.62) 



where 



m? m 



+ 



^Rj2{R„/2 + 1) 



y^Rj2{Rj2 + l) 



(7.63) 



and 



dv'^/Rj2{Rj2 + l) 



(7.64) 



where we defined ujq = K/m. To find solutions to equation Eq. (|7.62p we use the ansatz 



iP{b,r]) =e-y'^/^'^^{b,rj). 



(7.65) 



This leads to the equation for (j){b, rf) 



1 d'^cj) 
2m db"^ 



+ 



mco 



drj 



(7.66) 



The exact solution for late times is given by Eq. ()3.4p with initial conditions given by 
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Eq.dM]). Then Ea. dTIBS]) gives 

^ ^ g-OT/2m^(5^^) (7.67) 

where (p given in Eq. ()3.4p . 

As discussed in Chapter [3] an observer with a detector wih interpret the wavefunction 
of a given mode b at late times in terms of simple harmonic oscillator states at the final 
frequency u). 

In Fig. 17.81 we plot N versus v /Kg for various fixed values of u)Rs ■ We can see that the 
occupation number at any frequency increases as v/Rs decreases. Thus more particles 
are created as the shell reaches and crosses the horizon. However, the number of created 
particles does not diverge as R{y) — )• Rg. 

We then numerically evaluate the spectrum of mode occupation numbers at any finite 
time and show the results in Fig. 17.91 for several values of v/Rg. The first sign of non- 
thermality is the fact that the occupation number is non-divergent at w = 0, as opposed 
to the thermal Planck distribution in Eq. (j7.2ip . 

w7?5 = 0.15, 0.2, 0.5, 1 

N 




Figure 7.8.: Here we plot N versus v/ Rs for various fixed values of loRs- The curves are 
lower for higher values of u)Rs . 

In Fig. 17.101 we plot ln(l + versus ojRs for various values of v/Rg. As v/Rg 
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v/i?5 = -0.01,0.05,0.1 

N 




Figure 7.9.: Here we plot N versus ORg for various fixed values of v/Rg. The occupa- 
tion number at any frequency grows as the collapse progresses (i.e. v/Rg 
decreases) but in never diverges. 

decreases (as the shell is collapsing), the curves decrease. A thermal spectrum should 
gives us a straight line, however, we see that is not the case here. The best one can do 
is to fit the low frequency part of the spectrum and get the temperature in that regime. 
In our case we get T = {^.llRsY^. Unlike the case of Schwarzschild coordinates, where 
the spectrum becomes thermal in the whole frequency range, in Eddington-Finkelstein 
coordinates the spectrum never becomes thermal in the high frequency range. Another 
feature is apparent in Fig. 17.101 As the collapse progresses, the fluctuations in the 
spectrum become more violent. This is indicative of the shell approaching the actual 
singularity at ii = which is the region of strong gravitational fields. 

7.3. Discussion 

In this chapter we investigated the Hawking-like radiation produced during the time of 
gravitational collapse for both the asymptotic observer and the infalling observer. The 
occupation number of the radiation was then used to fit the temperature of the radiation 
as the shell approaches Rg. When considering Schwarzschild coordinates, in both cases 
the resulting analysis lead to the same conclusions: First, that the spectrum of the 
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y/Rs = -0.01, 0.001, 0.01, 0.05 

Ln(l + 1/N) 




Figure 7.10.: Here we plot ln(l + versus uRg for various fixed values of v/Rg. The 
curves are lower and display more fluctuations as v/Rs decreases. 

occupation display non-thermal characteristics during the time of collapse. This non- 
thermality is seen by a non-divergent occupation number when u) = Q and in oscillations 
about thermality. Second, the spectrum becomes more and more thermal as the domain 
wall approaches Rg, corresponding to large a) values. Finally, the spectrum becomes 
purely thermal when the domain wall reaches Rg- When consider Eddington-Finkelstein 
coordinates for the infalling observer, we find that the spectrum never becomes thermal 
in the high frequency range. 

In the case of the asymptotic observer, upon fitting the temperature, we find that 
the temperature of the radiation is on the order of the Hawking temperature. This 
value is not exactly the Hawking temperature for two reasons. First, when fitting the 
temperature we use a best fit approximation for the slope of /3. However, there is 
ambiguity for choosing the best fit approximation, thus the true slope of (3 may be 
different from the one chosen. Second, we are fitting the temperature numerically. There 
is always an inherent approximation used we using numerical methods, therefore our 
calculation is inherently ambiguous. 

In the case of the infalling observer, upon fitting the temperature in Schwarzschild 
coordinates, we find that when the shell reaches Rg the temperature of the radiation 
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becomes divergent. This would seem to imply that the local temperature measured 
by the observer is then infinite, meaning that the observer will burn up before he/she 
reaches Rg- However, this is not necessarily the case. It has been argued in Ref.[19j. 
where a simple 1 + 1 model was studied, that the local vacuum polarization will cancel 
out the divergent temperature energy density due to the radiation. Therefore, the true 
local value of the stress-energy tensor is small in the region R ^ Rg. A simple reason for 
this divergent temperature is that the Schwarzschild observer is actually an accelerated 
observer, so to truly investigate the local temperature one must consider a truly freely 
falling observer, i.e. the Eddington-Finkelstein coordinates. In this case, upon fitting 
the temperature, we find that the local temperature is in fact finite. 
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In 1972 Bekenstein argued that a black hole of mass M has an entropy proportional to 
its surface area, see Ref.[5D]. Further calculations by Gibbons and Hawking showed that 
the entropy of a black hole is always a constant, despite the type of metric which is used, 
see Ref. |21j. They showed that the expression for the entropy is given by 

Sbh = ^= ^Rl (8.1) 

where Ahor is the surface area of the event horizon and Rg is the Schwarzschild radius 
for a black hole which contains only mass. 

The typical method for calculating the entropy of the black hole is to first calculate 
the temperature of the black hole using the so-called Bogolyubov method. Here, one 
considers that the system starts in an asymptotically flat metric (typically Minkowski), 
then the system evolves to a new asymptotically flat metric (in the case of just mass, 
the typically final metric is that of Schwarzschild). One then matches the coefficients 
between the two asymptotically flat spaces at the beginning and end of the gravitational 
collapse. The mismatch of these two vacua gives the number of particles produced during 
the collapse. What happens in between is then beyond the scope of the Bogolyubov 
method, since the method is generally independent of time. Therefore the time-evolution 
of the thermodynamics properties of the collapse cannot be investigated in the context 
of the Bogoyubov method. 

Here we will investigate the time-evolution of a spherically symmetric infinitely shell 
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of collapsing matter in the context of the Functional Schrodinger formalism. Since 
the Functional Schrodinger formalism depends on the observer's degrees of freedom, 
one can introduce the "observer" time into the quantum mechanical processes, with 
the use of the Wheeler-de Witt equation, in the form of the Schrodinger equation, see 
Chapter [2j To study the case of gravitational collapse, one can then choose the classical 
Hamiltonian of the collapsing object, then employ the standard quantization condition. 
The wavefunctional is then dependent on the observer time chosen, hence one can view 
the quantum mechanical processes of a given system under any foliation of space-time 
that one chooses. The benefit of using the Functional Schrodinger formalism is that, in 
principle, one can solve the time-dependent wavefunctional equation exactly, as discussed 
in the previous chapters. Therefore the Functional Schrodinger formalism goes beyond 
the approximations of the Bogolyubov method, since the system is allowed to evolve 
over time, which allows one to investigate the intermediate regime during the collapse. 
Since the wavefunctional contains all the information of the system, one can, in principle, 
study the time evolution of the thermodynamical processes of the system. Of current 
interest is the time-evolution of the entropy of a collapsing gravitational object. We will 
do so from the view point of a stationary asymptotic observer, since this is the more 
relevant question. In this chapter we summarize the work originally done in Ref . |37] . 

8.1. Partition Function 

To study the entropy of the system, we will first develop the partition function for the 
system. In order to study the time-evolution of the entropy we shall employ the so-called 
Liouville-von Neumann approach, which was developed to study equilibrium and non- 
equilibrium quantum processes (see Ref. [26]). The Liouville-von Neumann approach is a 
canonical method which unifies the Liouville-von Neumann equation and the Functional 
Schrodinger equation. This approach utilizes the invariant operator approach developed 
by Lewis and Riesenfeld (see Ref. [2], Chapter [3] and Appendix [X]), which allows one 
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to exactly solve time-indepedent and time-dependent quantum systems. The Liouville- 
von Neumann approach has been employed for several different situations ranging from 
Condensed matter physics to Cosmology, see for example see Ref. |27j. The basic as- 
sumption of the Liouville-von Neuman approach is that non-equilibrium processes are 
consequences of underlying microscopic processes which are well described by quantum 
theory. The details about the collapse will depend on the particular foliation of space- 
time used to study the system. From the point of view of an infalling observer, in order 
to calculate the backreaction and local effect around the event horizon it is important to 
choose a state that is non-singular at the horizon. In this region, the vacuum of choice 
is the Unruh vacuum (see Refs.|28t I29|). However, discussed above, we are interested in 
the view point of the asymptotic observer. 

Using the Liouville-von Neumann approach, and following the procedure used in 
Ref. [27], we can write the partition function as 



and /? is a free parameter. Here we note that Eq. (j8.3p is just the Heisenberg equation of 
motion for the operator I, see Ref. [22], where the total time derivative of the operator 
is zero. In the case that the total derivative is equal to zero in the Heisenberg, this case 
is known as the Liouville-von Neumann equation, see Ref. [2]. 
From Ref. [3], we can write the invariant operator / as 



Z = Ti- e~^' 




where I is any operator which satisfies the equation 
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Here we note that the invariant operator / is time dependent since p is time dependent 
(see Eq. ()3.5p ). Using Eq. ()8.4p we can therefore write the partition function, Eq. ()8.2p . as 



Z = Ti exp 



'^2 



- + (vTb/O - mprjby 



.5) 



In this form we can see that the partition function is time dependent since the invariant 
operator / is time dependent by virtue of Eq. ()8.4p . 

We note that we can rewrite the invariant operator in a more suggestive manner by 
writing Eq. (|8.4p as 



I = ^ 



1 

V2 



nit) + - 



p) 



i (irbp - mprjb) 



-\ + i (vTbp - mp^6) 



(8.6) 



where 



n{t) = a\t)a{t) 



and 



a{t) = 



V2 



1/4 



+ i {TTbP - mpnb) 



(8.8) 



Here n{t) is the time dependent number of states. Hence, the invariant operator / takes 
on the form of a time-dependent harmonic oscillator Hamiltonian, where the number 
operator is time dependent. 

For a physical meaning of the partition function, we need to act the invariant operator 
on a quantum state. In the Heisenberg picture, the quantum states span a particular 
Hilbert space. A convienient basis in this Hilbert space is the so-called Fock space 
representation, see Ref . j23] . This basis is an eigenstate of the Number operator, Eq. ([87 
Thus at a particular time t, one has in the Fock space representation 



n{t)\n, t) = n\n, t). 



(8.9) 



69 



CHAPTERS. ENTROPY 



Thus, in this space we can then write the partition function as 

Trexp 
1 



-/3ujo [n + ^ 



2 sinh 



(8.10) 



At first glance, one would be tempted to say that the partition function in Eq. (j8.10p is not 
time-dependent since the partition function now only depends on the initial frequency of 
the induced scalar field. However, recall that /3 is free parameter which we can choose. 
Here we discuss our choice in the free parameter /3. 

In Refs.[30t [32] one can define the occupation number for a frequency O, Eq. ()3.2ip . 
Then by fitting the number of particles created as the usual Planck distribution Eq. ()7.2ip . 
one can then in principle fit the temperature of the radiation. Here, we then choose to 
define f3 as 

/, = ^Ml±iM. ,8.11) 

OOJ 

This implies that all of the time dependence of the system is encoded into the tempera- 
ture of the system. 

Therefore we can see that Eq. ()8.10p is just the standard entropy for a time-independent 
harmonic oscillator, however, the temperature here is time-dependent. Thus we recover 
the time-dependence of the partition function. Since the partition function is time- 
dependent, therefore the entropy is also time-dependent. 

8.2. Entropy 

In terms of the partition function, the thermodynamic definition of entropy is given by, 
see for example Ref. |24j . 

S = lnZ-(3^. (8.12) 
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Using Eq. (j8.10p . we can then write the entropy of the system as 




(8.13) 



1 _ g-Z^'^o ■ 



Therefore, this is again just the entropy of the usual time-independent harmonic oscil- 
lator. Prom Eq. ()8.1ip it follows that the temperature is time-dependent. 

To be able to calculate the entropy of the domain wall we will consider the entropy of 
the entire system, i.e. the domain wall and radiation, and the radiation alone. We will 
assume that the total entropy is a linear equation in the entropy of the domain wall and 
the entropy of the radiation. Thus we will write the total entropy as 



where the subscripts SR stands for domain wall and radiation, S for just domain wall 
and R radiation only, respectively. Then by subtracting these two quantities one can 
then determine the entropy of the domain wall 



In Chapter \7\ we considered the wavefunction and occupation number of the radiation 
only system. To proceed further, we must now consider the wavefunction and occupation 
number for the entire system, SR. 

8.2.1. Entire System 

To find the wavefunction and occupation number for the entire system, we first note 
that from Eq. (|5.5p we can approximate the Hamiltonian of the domain wall as 



SsR = Ss + Sr 




Ss = SsR - Sr. 




BUr. 
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Then using Eq. (|8.16p and Eq. (j7.13p we can write the Hamiltonian of the entire system 
as 

H = H^aii + Hk = -BUR + B^ + ^b^ (8.17) 

2m 2 

where Hji is given in Chapter [5] and Uf, is given by 

n, = -i^ (8.18) 

The wavefunction for the entire system is then a function of 5, R, and t, which we can 
write as 

= '^{b,R,t). (8.19) 

Substituting Eq. (j8.17p into Eq. (j2.12p . we can then write the Functional Schrodinger 
equation as 

To solve Eq. (|8.20p we wih use the semiclassical case, i.e. we will use the classical back- 
ground for the collapsing shell. Since the distance of the shell only depends on the time, 
see Eq. (|5.13p . we can then write 



dt B d'^^ o 
\ /)^\& = i 

dR dt 2m 562 ^2 dt' 



Hence, we are eliminating the R dependence from Eq. ()8.19p . so ^(b,R,t) — )■ ^{b,t). 
Rewriting gives 

B d"^^ K ^ d^ f dt\ , , 

+ —b^q, = i [l- B—] . (8.21) 



2m db"^ 2 dt \ dR J 

Making use of Eq. ()5.13p . i.e. dt/dR = —B, this becomes 



B d"^^ d^ 

+ —b^qr = 2i—. (8.22) 



2m db'^ 2 dt 
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We now rewrite Eq. ()8.22p in the standard form 



1 9^ m 2/~\,2 
l^d^^j"" ^^^^ 



^(M) = ^^ (8.23) 



where 



and 



fl = l tdt'(l-^] (8.24) 



2 ./o V R 



ml — Rs/R 1 — Rs/R 

Here we have chosen to set f/(t = 0) = 0. 

The solution to Eq. (|8.23p is given by Eq. (|3.4p , as discussed in Chapter [71 We can then 
find the occupation number N for the entire system, Eq. (j3.2ip . 

Here we will make some quick comments regarding the occupation number. We can 
see that from Eqs. (|8.23p and ()7.14p . the Schrodinger equations for the entire system and 
radiation only are of the same form. Hence one would expect that there is no difference 
between the occupation number for the entire system and the radiation only. However, 
the time parameters fj and i], given in Eqs. ()8.24p and ()7.15p . are different. Hence the 
occupation numbers of the two systems will evolve differently, which leads to different 
temperatures in each of the two systems. Therefore, the entropy of each system will be 
different. 



8.3. Analysis 

First we consider the entropy of the entire system. In Figure \8A\ we plot the entropy of 
the entire system as a function of dimensionless time t/Rg. Figure 18.11 shows that the 
system starts with an initial entropy of zero. This is expected since initially there is only 
one degree of freedom, meaning that S = ln(l) = 0. Here we have normalized the initial 
entropy of the shell to be zero. To justify this normalization, consider a solar mass black 
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hole. Under the usual Bekenstein-Hawking entropy, the order of magnitude estimate 
of the entropy of a solar mass black hole is Sbh ~ lO"^^. Now consider that the shell 
is actually made up of protons. The initial entropy of the shell then is approximately 
Ss,o ~ 10^'^. Comparing the entropy of the final black hole versus the initial entropy of 
the shell, the entropy of the final black hole is much much greater than that of the initial 
entropy of the shell, thus the initial entropy of the shell only contributes a negligible 
amount of entropy to the entropy of the final black hole. Thus our normalization of the 
initial entropy of the shell to zero is justified. As t/Rs increases, initially the entropy 
increases rapidly, then settles down to increase approximately linearly. Due to the linear 
increase, we see that as t/Rg goes to infinity, the entropy will then diverge. This is again 
expected since as the asymptotic time goes to infinity, the number of particles that are 
produced diverges (see Ref . |30j ) . This is a consequence of the fact that we keep the 
background fixed (i.e. Rg is a constant). In reality, Rg should decrease over time since 
the radiation is taking away mass and energy from the system. Therefore as t/Rg goes 
to infinity, the entropy of the entire system as measured by the asymptotic observer 
diverges as R ^ Rg. 




5 10 15 20 25 30 

Figure 8.1.: We plot the entropy of the entire system as a function of asymptotic observer 
time t. 

This is consistent with the results found in Refs.[33]. Here the authors consider the 
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time-dependent non-equilibrium evolution of a black hole as well as the incorporation of 
the given off radiation. Here one can see that the entropy of the system diverges as the 
time goes to infinity. 

The results of Figure 18.11 are consistent with the generalized second law of black hole 
thermodynamics. The generalized second law states that, see for example Ref . |38j and 
references there in 

5{Sout + A/4) > (8.26) 

where here, Sout = Sr, A/A = Ss and Sout + A/4 = Ssr, respectively. Eq. (j8.26p 
simply states that the total entropy of the system must constantly be increasing as in 
agreement with thermodynamics entropy Ref. [2^. As stated above, a realistic model 
for gravitational collapse will have that the Schwarzschild radius Kg will decrease over 
time, since the domain wall is losing mass. Eq. ()8.26p allows for this result as long as the 
entropy increase of the radiation compensates for the loss in entropy of the collapsing 
domain wall. 

Now we consider the radiation only. Considering just the particles which are created, 
i.e. the radiation, during the collapse, we can then plot the entropy as a function or the 
rescaled asymptotic time t/Rg, see Figure Figure [52] shows initially the entropy of 
the system is zero. Again, this is expected since initially the domain wall is in vacuum, 
meaning that there are no particles produced. Therefore the only degree of freedom is 
that of the domain wall, this then gives that the initial entropy must be zero. As the 
asymptotic observer time increases, initially there is rapid increase in the entropy, but 
again, the entropy then increases linearly as the asymptotic observer time increases. As 
in the case of the entire system, as the time measured by the asymptotic observer goes 
to infinity, the entropy of the particles created during the time of collapse diverges. This 
is expected since the number of particles which are created during the time of collapse 
diverges as ii — )• Rs, hence as the domain wall approaches the horizon the number of 
particles created during the collapse diverges. This result again is in agreement with the 
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generalized second law of black hole thermodynamics, Eq. (j8.26p . 

S(t) 

30 r 



25 
20 
15 
10 
5 



10 15 20 25 30 



t/Rs 



Figure 8.2.: We plot the entropy of the particles created during the collapse as a function 
of asymptotic time t. 



In Figure 18.31 we plot the entropy as a function of the rescaled asymptotic observer 
time t/Rs of both the entire system and the particles created during the time of collapse. 
Figure [8^3] shows that except for the initial increase in the entropy, for later asymptotic 
observer time, the slopes of the entropy versus time are approximately equal. Therefore, 
one can expect that the entropy of the domain wall is approximately constant for late 
times. 




5 10 15 20 25 30 



t/Rs 



Figure 8.3.: We plot the entropy as a function of asymptotic observer time t for both 
the entire system and the particles created during the time of collapse. 
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As stated earlier, what is of interest is the entropy of the collapsing domain wall, 
since this will collapse to form a black hole. To find the entropy of the domain wall, we 
can take the entropy of the entire system and subtract off the entropy of the particles 
produced (since these are the only relevant objects which contribute to the entropy), 
see Eq. ()9.3p . The result is then given in Figure 18.41 Figure 18.41 shows that initially 
the entropy of the domain wall is zero. As stated above, this is expected since initially 
there is only one degree of freedom. As asymptotic time increases, the entropy of the 
domain wall rapidly increases. However, for late times, the entropy of the domain wall 
goes to a constant. As stated above, this is expected since the late time entropies for 
entire system and for the particles created during collapse are approximately parallel. 
However, as discussed earlier, one would expect that in a realistic model the entropy 
of the domain wall should in fact decrease over time since Rg is decreasing because the 
domain wall is losing mass. The entropy here, however, is constant since we are assuming 
that the mass is approximately the Hamiltonian of the system, which is a constant of 
motion, see Chapters [5] and HI This means that since we are holding the mass of the 
domain wall constant, we need to keep adding energy to the system to counter act the 
loss of mass from the Hawking radiation. Therefore one can expect that the entropy of 
the domain wall must be a constant for late times. 

In reality, radiation takes mass away from the system, so the entropy of the domain 
wall will go to zero as Rg goes to zero. This means that after the black hole disappears, 
all the entropy will go into the entropy of the radiation, which is in agreement with the 
generalized second law of black hole thermodynamics. 

From Figure 18. 4|, we see that our numerical value for the late time entropy of the 
domain wall is 

S ^ 0.7i?^ 

Comparing with Eq. (j8.ip . we can view this discrepancy as a shift in the Schwarzschild 
radius Rs- In order to get the theoretical value for the entropy, Eq. (|8.ip . we see that we 
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Entropy of Shell 




Figure 8.4.: We plot the entropy of the shell as a function of asymptotic observer time t. 

would require Rs — )■ 2.11i?s- This is an understandable numerical error, which implies 
that our numerical solution is of the same order as the Hawking-Bekenstein entropy. 

Another interesting thing to note is that Figure 18.41 tells us that change in entropy 
occurs for early times, then gets frozen as time increases. From the plot we see that the 
change in entropy occurs during the time range < t/Rg < 7.5. At first sight this seems 
to be an arbitrary value for the entropy of the domain wall to stop increasing. However, 
from Eq. (j5.14p one can see that this time is not an arbitrary value. 

To see this, let us first consider Eq. ()5.14p and make the requirement that Rq = nRg, 
where n is some integer. Then we can write Eq. ()5.14p as 

i?(t) = i?, (l + (n-l)e-*/^=). 

For illustration purposes let's restrict the value of n to be n < 10, which is a restriction 
that the domain wall starts off at a position ten times it's Schwarzschild radius. In 
Figure [831 we plot R/ Rs versus t/Rg for various values of n. For each value of n chosen, 
we see that the value R/Rg ~ 1 occurs for t/Rg ~ 7. In the case of n = 10, we see 
that R = l.OObRg, while the value is less than that for smaller values of n. Hence, the 
time t/Rg = 7.5 seems to be a universal time when the domain wall is almost to the 
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Schwarzschild radius. From Eq. (j5.13p we see that by this time we have 

R=-B^O. 

Hence in this time hmit, the velocity of the domain wall is approximately zero, meaning 
that as far as the asymptotic observer is concerned the domain wall has stopped moving 
and there are no more dynamics. This can be seen in Figure 18.61 where we plot the 
corresponding velocities for the same values of n. Figure 18.61 also shows that the time 
t/Rg = 7.5 corresponds to a universal time of when the different velocities go approxi- 
mately to zero. Recall from Chapter O that it takes an infinite amount of time for the 
domain wall to reach the horizon, so from t/Rg = 7.5 to infinity the entropy is constant 
since all the dynamics are essentially done and the shell is approximately stationary 
for the observer. Hence the volume of the spherically symmetric domain wall becomes 
essentially constant by the time t/Rg = 7.5. 

n=2, 5, 10 

R/Rs 




Figure 8.5.: We plot R/Rg versus t/Rg for various values of n. Here the blue curve 
corresponds to n = 2, the green curve corresponds to n = 5 and the red 
curve corresponds to n = 10. 

Second, we can show that the entire system and the induced radiation come into ther- 
mal equilibrium at this time. In Figure [8.1 II we plot /3 versus t/Rg for the entire system 
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(continuous curve) and the induced radiation (dashed curve) . Figure 18.111 shows that for 
the time t/Rs ~ 7.5 the values of the two /3's become approximately equal, meaning that 
the entire system and the induced radiation are now at the same temperature. There- 
fore the system is now in thermal equilibrium, meaning that there is no more change in 
entropy of the domain wall as t/Rs increases. Further more, the fluctuations (departure 
from thermality) in /? become very small at this time, as discussed in Refs.|30l I32j. 

n=2, 5, 10 



(dR/dt)//?5 




Figure 8.6.: We plot R/Rg versus t/Rg for various values of n. Here the blue curve 
corresponds to n = 2, the green curve corresponds to n = 5 and the red 
curve corresponds to n = 10. 

Finally we can evaluate the the chemical potential for both the entire system and for 
the induced radiation. From definition we can write the chemical potential as 

.-f. (8.27) 

In Figure 18.71 we plot the chemical potential for both the entire system and for the 
induced radiation. We can see that as t/Rg increases the chemical potential of the 
entire system and the induced radiation goes to zero. This means that the dispersion of 
particles goes to zero and the system goes into equilibrium. 

During the dynamical process, the entropy increases almost linearly. If one applies a 
best-fit line, we see that the entropy oscillates about the best line. These oscillations may 
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Figure 8.7.: We plot /i versus t/Rg. The solid line corresponds to the entire system while 
the dashed line corresponds to the induced radiation only. Here we see that 
as t/Rs increases, the chemical potential for each goes to zero. 

be attributed to several different circumstances. First, the oscillations may be caused 
by the non-thermal property of the radiation (see Ref . |30j ) . Secondly, these oscillations 
may be a manifestation of the error associated with the numerical calculations. Lastly, 
the oscillations may be an artifact of expanding the calculations beyond the region of 
validity, since we are using the near horizon approximation. Hence for values large 
compared to Rg, we cannot completely trust our result. 

8.4. Discussion 

Here we have shown that the entropy of the collapsing domain wall and the entropy of 
the radiation given off during the time of collapse are in agreement with the generalized 
second law of black hole thermodynamics. The results of Figure 18.11 are clearly in 
agreement with Eq. ()8.26p . The results of Figure 18.41 are in agreement with the results of 
Hawking and Gibbons, Eq. (j8.ip . that the entropy of the black hole is in fact finite and 
proportional to the area of the event horizon. 

Note, here we do not discuss or explain the origin of Eq. ()8.ip . we merely verify that our 
model gives the correct result. The origin of Eq. (|8.ip is still not understood, however. 
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many attempts have been made to make sense of this result (see for example Refs. [3^l35^ 
I36j). However, the answer to this question may lie in understanding the entanglement 
nature between the particles inside and outside of the event horizon, see for example 
Refs.lMlllQiSI]. 
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In this section we make some general comments on how one can include the effect of back 
reaction for an infalling observer during gravitational collapse. In this section we do not 
completely solve the equations of motion for the included back reaction, we merely set 
up the situation and make some comments about it. 

To incorporate back reaction into gravitational collapse, one must consider the entire 
Hamiltonian, as in Chapter[8l as well as the interaction Hamiltonian between the domain 
wall and the induced radiation. 

Thus the total Hamiltonian is given by 



H =Hwall + Hjiad + Hint 



modes 



2m 2B 



+ T^,S^'' (9.1) 



where Hjnt = T^^S^^^ is the interaction Hamiltonian and T^y and 5^^^ are the energy- 
momentum tensors for the radiation and the domain wall, respectively, which are given 

by 



r, 



1 



27r2 dt drr^ 



(9.2) 
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and 



(9.3) 



From the expansion of the scalar field in Eq. ()3.8p . we can see that the stress-energy 
tensor for the scalar field takes on the form 



Too 


Toi 








Tio 


Tn 














T22 














^33 



(9.4) 



While from Eq. (|9.3|) we see that the stress-energy tensor for the domain wall takes the 
form, 



goo 











5" 























(9.5) 



Hence we see that the interaction Hamiltonian doesn't contain any off-diagonal terms. 



9.1. Stress- Energy Tensor 

Here we develop the stress-energy tensor for the radiation and the domain wall, respec- 
tively. First we will discuss the stress-energy tensor for the induced radiation. Second 
we will discuss the stress-energy tensor for the domain wall. 
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9.1.1. Radiation Stress-Energy Tensor 

Here we examine the stress-energy tensor for the radiation. From the discussion in 
Chapter [7] we can write the stress-energy tensor as 



=47r / (It 



drr' 



B {dr^f 



+ / drr 
Jo 



Rir) ^JB + Rl l-Rs/r 

2^ 



J Rir ) B V r J 



drr'^d^^dy^ 



(9.6) 



Now using the metric we can write the individual terms, which are given as 



Too =47r / dT 



B 



drr^ (1 ^ I — , 

r J .JlTM 

R{r) 



drr , 



■r Jo 



drr 1 



R 



Jr{t) B \ r J ^ Jo 



(9.7) 



Toi = -in / dr / drr^dr<^dr<^, 



-10 



4iT dr drr'^dr^dr^, 



(9.8) 
(9.9) 



Til =4vr / dr 



1 
2 

drr^ 



R{t) 



drr' 



B 



+ 



R{r) 

CO 

R{r) 



^/BTM{l-Rs/rf 



+ / drr'^^/l + Rl- 
Jo ^ 



Rs/r 



drr' VB + R'r ^g^^^2 
B 



drr'^{dr^f 



(9.10) 
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T22 =47r J dT 

r.R(r) 



1 



R{r) 

drr — , 

^/TTM 



00 

drr^ — , — 

R[r) y/B + R^l- Rs/r 



Jo Jr{t) ^ V r J 



and 



^33 =2vr y dT 

Rir) 



R{r) 



drr 







B 



drr"^ — , — — — 

R[r) y/B + R^l- Rs/r 



+ / drr^y/T+W,{dr<^f + r drr^^^?^^ fl " — ) (^r^)') • (9.12) 
/o Jr{t) 



B 



Here note that from Eqs. (f9TT]l and (f9l^ show that T33 = {2/3)T22. 

For a full analysis of the stress-energy tensor we will look in the near the horizon 
limit. Ideally we would like to extend this analysis to the near singularity limit as well. 
However, we are working in Schwarzschild coordinates, which we cannot extend to the 
near singularity limit due to the fact that the observer is being constantly accelerated 
(see Chapter [5]) . We will then exam the behavior of the stress-energy tensor near the 
horizon, i.e. in the region i? ~ i?^. 

Of interest is the behavior of the stress-energy tensor near the horizon. To investigate 
the effect of the radiation we will change the limit of integration from R{t) to Rs, 
allowing us to find the dominate terms in this regime. Prom Eq. (|9.7p we see that in this 
limit and with the expansion in modes we have 



'-00 



47r / dT 



Rs 



drr^ f 1 - ^ 







+ 



Rs 



dT 



B 



r 



drr^{dr<^f 



1 1 . I . 1 yjB + Rl ~ 1 . . 

n / =0'k-t^kk'ak' — 5 akCkk'O-k' — -^ak^-^ kk' ^k' 

2 ^l + R^ IB I 
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from Eg. ([97 



drdk^kk'ak' 



from Eq.dM 



from Eg. (19111 



T 



11 



drr 



1 



poo 

+ / drr 



2^/B + Rl,^ 2\ /"°°^_2/o^^2 



-(3.<l>)' 



drr^((9^^>)" 
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and from Eq. (|9.1ip we have 
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where the matrices are defined by 



2 I - 



Akk' =47r / drr' U " — fkfk', (9.13) 



r 

OO / 7~> \ 2 



OO 

2 fl 



Dfcfc, =87r / drr'f'Jk,, (9.15) 







OO 



Efcfc, =47r / drr^hfk,^ (9.16) 







47r / cirr2 ( 1 - i::! ) 7;,/^,, (9.17) 







r 



OO 



Gfcfc' =47r / drrVfc/fc', (9.18) 



OO 



Hfcfe, =87r / drr'f'J'j,,, (9.19) 
) 



Jfcfc, =47r / drr^fkfk', (9.20) 

JO 

Kfcfc' =47r^"dr/ (^1 - ^) /(/^,. (9.21) 

To further investigate the problem, for a moment let us assume that the basis functions 
are planewaves. This is a vahd approximation in the asymptotic regime, however, this 
will give us some insight into the problem here. For the basis functions as plane waves 
we have 

£ ikr 
Jk = e , 

however since we are requiring real basis functions then we will take the real part of this. 
Therefore we can write 



(/•OO \ / POO 

J drr^f'jk'] = ReU drr'^e'^^+^'> 
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Performing the integral over r we then have 



Re 



\ 

drr^f'Jk'j = Sik + k') 



which is finite. Hence the other terms in Eqs. (|9.7p . (|9.10p . (|9.1ip and ()9.12p are dominate 
due to the divergences of these terms. Therefore we can ignore these extra terms, so we 
then have 



Too 

Toi 
Tio 

Til 
T22 



dr 



-j^ ^ 'Q _|_ ^2 

+ i?^dfc(A"^)fcfc'dfc/ - — -akCkk'dk' 

I Id 



drak^kk'ak' 



drakBf^l/dk', 



B 



dr 



-y/l + R'^dk{J~^)kk'dk' + -akKkk'CLk' 



which is of the same form as Hrad- We can see that the matrices in Eq. ()9.2ip are just 
multiples of the matrices in Eq. (|7.3ip . Therefore we can see that the eigenvalues of 
Eq. ()9.2ip are multiples of those of Eq. (|7.3ip . hence we can simultaneously diagonalize 
the matrices as we did in Chapter [71 Finally we can write 



-too 

Toi 
Tio 

Tu 
T22 



dr 



1 VTT^ 5^ _ nI£y^BTM,2 
'2 nm dh"^ 2 B 



drdk^kk'^k', 



draki^f^^idk', 



dr 



dr 



2 m 562 2n 5 



2 m 962 2 B 



(9.22) 

(9.23) 
(9.24) 

(9.25) 
(9.26) 
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where as in Chapter[71 m and K are eigenvalues, b are the eigenmodes and n is a constant 
multiple. Here we note that Too, Tn and T22 have the structure of a Harmonic oscillator. 

We now calculate the expectation value of the stress-energy tensor. To do this we 
consider 

{T,,) = {0\T,,\0) (9.27) 

where |0) is the vacuum state. Since the components of the stress-energy tensor have 
the structure of a harmonic oscillator we take that the vacuum state is the ground state 
of the harmonic oscillator (see Chapter [7|) . The ground state of the harmonic oscillator 
is given by Eq. ()7.18p . thus the expectation value is 

(T^^) =Jdb ,-rn.„6V2j.^^^-™a.o6V2. (g ^S) 

From the structure of Eqs. (j9.22p - (j9.26p . we can see that there is a kinetic term and 
a potential term which all have the same dependence on the eigenmode b. So, using 
Eq. ()7.18p we can write 



Kinetic Term 



362 



4 ' 

and 

Potential Term = J db (^-^^^^'' 1%"^ ^r'^^^^'' 



where we used the fact that there are an infinite number of eigenmodes b (hence the 
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integrals are from zero to infinity). The individual components are then, 



(Tn) 
{T22) 



dr 



dr 



dr 



2 nm \ A ) 2 B \A{muQfl'^ ) 

I n^JlTM _ K ^B + Rl / ^ 

2 m \ A ) 2n B V^C^nwo)^/^ 



'2 m V 4 y 2 B \A{mojQYl^) 



(9.29) 
(9.30) 
(9.31) 



Investigating Eqs. ()9.29p - ()9.3ip we can see that as R ^ Rg the potential term diverges 
such as in the Hamiltonian of the induced radiation. Therefore we can conclude that 
the components of the stress-energy tensor are divergent, however, not due to the usual 
reasons. Typically this divergence of the stress-energy tensor is associated with the 
divergence in the frequency u (see for example Refs.[28l \29\ 123]). To get around this 
divergence, one usually either applies a cut-off for the allowed frequency or applies a 
renormalization technique that makes the stress-energy tensor finite. Here we can see 
that this process is not needed since the divergence is not due to the frequency (since 
we never specify the basis functions), however the divergence is due to the metric itself. 

The divergence in Eqs. ()9.29|) - ()9.3ip in the regime ~ i?s is due to the B term in 
the potential term. As stated earlier, Chapters [5] and [71 this is due to the fact that we 
are using Schwarzschild coordinates. The Schwarzschild observer is in an accelerated 
reference frame, which causes the divergence. Therefore, as we saw in Chapter [71 to 
study the question of backreaction the more appropriate observer to use would be a 
truly free-falling observer such as an Eddington-Finkelstein observer. However, we will 
not investigate such an observer here. 
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9.1.2. Stress-Energy of the domain wall 

From Eq. (|9.3p we can write the determinant of the induced metric as 



-j = R^sm9jB-^ 



where we used Eq. (j4.7p . The stress-energy tensor for the domain can then be written as 

S^"' = -a-i"^ (9.33) 

where 7'''' is again the induced metric on the surface of the domain waU. This is expected 
from Eq. (|B.10p in the case of the domain wah [a = rf). 

9.2. Quantum Hamiltonian 

Here we wish to find an appropriate way to take into account the fact that the mass of 
the domain wah is changing, due to the fact that the radiation is taking mass away from 
the system. To do this we wih follow a technique that was first introduced in Ref . |43j . 

Eq. ()4.12p tells us that the mass of the domain wall is approximately the Hamiltonian, 
therefore we can write the Schwarzschild radius as 

Rs ^ 2GHwaii- (9.34) 

The factor B in the Hamiltonian for the radiation contains the energy of the wall via 
the Schwarzschild radius. So we then have 



B = l- (9.35) 
it 
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In the near horizon Hmit Eq. ()4.12p can be written as 



M = 4^c7fi2[^i + i?2 _ 27rGaRs] = H^aii- (9.36) 

Assuming that the velocity at the horizon is smaU and dropping the constant terms we 
can write this as 

H^aii = 2^aRlRl. (9.37) 
Using Eg . (19.341) we can then rewrite this as 



(9.39) 



where in the second hne we used Eq. ()5.17p . 

The total Hamiltonian in terms of a single mode then becomes 



J L T S^'^ (9 401 

2 {{IQiraG^y/m - 2G(n2j)V3) (ul)2/3 ' ^^'^^^ 

Here we note an unusual property of Eq. ()9.40p . the appearance of the fractional deriva- 
tives. In general, fractional derivatives are non-local, that is, one cannot say that the 
fractional derivative at a point x of a function / depends only on the graph of / very 
near x, see for example Ref.[Mj. Therefore it is expected that the theory of fractional 
derivatives involves some sort of boundary conditions, involving information further out. 
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The most general definition of tlie fractional derivative is 



aJ^t 



dxi 1 



Re{q) > 
1, Re{a) = 

Jl^idx)-", Re{a) < 0. 



Here the first case is defined as 



dxi' 



(k-a) 



■X 



k—q 



r(fc + i) 

T{k-q + iy 



while the third case is defined as 



(9.41) 



(9.42) 



(dx)l^l/(x 



{x-tt\-^f{t)dt. 



(9.43) 



Hence Eq. (|9.4U|) is an differential- integral equation. 

The study of the behavior backreaction is therefore very complicated. However, the 
interesting thing to point out here is that, similarly to the investigation of the quantum 
mechanical effects studied in Chapter El the presence of the non-locality again emerges. 
However, in Chapter [6] the non-locality was only present when investigating the near 
classical singularity regime. Here, the non-local effect is even present in the near horizon 
regime. 



9.3. Discussion 

In this chapter we investigated the stress-energy tensor for the radiation given off during 
collapse as well as investigated a way to include the loss of mass during this collapse. 
We found some interesting properties of these two quantities. 

First for the stress-energy tensor for the radiation, we found that the expectation 
value for the stress-energy tensor is in fact infinite. However, this is not due to the usual 
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difficulties. Generally when one investigates the stress-energy tensor, the infinities arise 
from the basis function. Traditionally one assumes a plane-wave basis function for the 
radiation, and the divergence is therefore due to the frequency of the basis function. To 
avoid these infinities, one usually institutes a cut-off frequency. Here, we do not have 
this problem. This is due to the fact that we never actually specify our basis functions, 
hence we do not have the problem of infinities in the basis function. The divergence in 
this case is due to the presence of the in the potential term. As — t- Rs, i? — > 
which causes the divergence. As stated in Chapters [5] and [7] this is due to the fact that 
we are using Schwarzschild coordinates, where the observer is being accelerated. 

To include the loss of mass into the Hamiltonian of the system, we used the technique 
originally developed in Ref . |43j . Here one uses the approximation that the Hamiltonian 
of the domain wall is approximately the mass of the domain wall. Therefore one can 
replace the mass in the Schwarzschild radius by the Hamiltonian of the wall. Using this, 
we can then rewrite the total Hamiltonian as in Eq. ()9.40p . The interesting thing here 
is that the Hamiltonian is now in terms of fraction, not whole or partial, derivatives. 
By definition, fractional derivatives are not strictly local quantities and will either give 
a differential or integral equation depending on the sign of the fractional derivative. As 
in Chapter [6l we recover the non-locality of the quantum effects during gravitational 
collapse. Unlike in Chapter El these effects are now manifest even near the horizon. 
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In this thesis we have investigated quantum mechanical effects of gravitational collapse by 
utilizing the time-dependent nature of the Functional Schrodinger formalism. As stressed 
throughout this thesis, the Functional Schrodinger formalism allows us to investigate 
the intermediate regimes that the standard methods cannot. Therefore we can obtain a 
better understanding of what is happening during the evolution of the collapse, at least 
in the context of the Functional Schrodinger formalism. As we have seen in the previous 
chapters, the effects in this intermediate regime are robust and give good insight into 
the process of collapse. 

This thesis is not meant to be an exhaustive list of the different types of gravitational 
collapse. Here we solely concentrated on a massive domain wall, while ignoring all other 
observable quantities (such as charge and angular momentum). However, one can "eas- 
ily" incorporate these observable into the system as well. For example, one can repeat 
the steps above for the case of a massive-charged domain wall (i.e. Reissner-Nordstrom) . 
This has been done for the classical and quantum solutions in Ref . |45j and for the semi- 
classical radiation in Ref. |46j . The analysis can also be repeated for different topologies, 
other than spherically symmetric domain walls, as well as for different asymptotic space- 
times (such as de Sitter or anti-de Sitter). In Ref. |47j the classical and quantum solutions 
are studied for a (3-|- l)-dimensional BTZ black string in AdS space. It is well know that 
a (3 -|- l)-dimensional BTZ black string has the topology of a cylinder and is asymptotic 
to AdS space-time, due to the negative cosmological constant. 
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It is also important to note here that the Functional Schrodinger formalism is not 
restricted to gravitational collapse. One could also apply the formalism to expanding 
systems as well, which are essentially collapsing systems in reverse. In this case, one can 
investigate an expanding de Sitter or anti-de Sitter universe and consider the radiation 
and entropy during the evolution of expansion. For subsequent work on the radiation 
given off during expansion, see "Time dependent fluctuations and particle production 
in cosmological de Sitter and anti-de Sitter spaces," by E. Greenwood, D. Dai and D. 
Stojkovic (submitted for publication in Phys. Rev. D). In the case of de Sitter expan- 
sion, which is represented by the Freedman-Robertson- Walker metric, the horizon is the 
largest comoving distance which light emitted now can reach the observer at any time in 
the future. It is expected that that de Sitter space can produce thermal radiation as well 
(for some counter arguments see Refs.flH} I49j). In the case of anti-de Sitter expansion, 
unlike de Sitter expansion, the space-time does not contain an event horizon. Therefore, 
one would not expect thermal radiation with a constant temperature. However, due to 
the time-dependent metric, particle production is still expected. Here it is expected that 
after a short time of expansion, the universe starts recollapsing and ends up forming a 
black hole, see for example Refs.[50l I51j. The Functional Schrodinger formalism can be 
applied to these situations as well to help shed light on these questions. 
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In this section we breifly review the invariant operator method developed by Lewis and 
Reisenfeld in Ref.[2j as a solution to the time-dependent Schrodinger Equation. 

Consider a system whose Hamiltonian operator H{t) is an explicit function of time, 
and assume the existence of another explicitly time-dependent non-trivial Hermitian 
operator I{t), which is invariant. To say that I{t) is invariant means that I{t) satisfies 
the Liouville-von Neumann equation 



Here we will consider the analysis for a state vector ^tp), however, in general this also 
works for a wavefunction since 'iIj{x) = {x\ip). We can then write the Schrodinger equa- 
tion as 



By operating with the left-hand side of Eq. (IA.l| ) on the state vector and using the 




(A.l) 



and since I{t) is Hermitian we have 



/t = /. 



(A.2) 




(A.3) 
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Schrodinger equation, we obtain the relation 



.d 



{I\^P)) = H {I\^)) 



(A.4) 




which imphes that the action of the invariant operator on a Schrodinger state vector pro- 
duces another solution of the Schrodinger equation. In general, this result is valid for any 
invariant, even if the invariant involves the operation of time differentiation. However, 
for our purposes, we shall consider invariants which do not involve time differentiation. 
This choice allows one to derive simple and explicit rules for choosing the phases of the 
eigenstates of I{t) such that these states themselves satisfy the Schrodginer equation. 

Assume that the invariant is one of a complete set of commuting observables, so that 
there is a complete set of eigenstates of /. Denote the eigenvalues of / by A, and the 
orthonormal eigenstates associated with a given A by |A,k), where n represents all of 
the quantum numbers other than A that are necessary for specifying the eigenstates: 



Since the invariant is Hermitian, the eigenvalues A are real. They are also time-independent 
as we shall now see. By differentiating Eq. (|A.5p with respect to time, we obtain 




(A.5) 




(A.6) 




(A.7) 



Using Eq. ()A.ip we can write 
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The scalar product of Eq. (jA.8p with a state I A', k') is 



which then imphes 



.81. 

i{X' , k'\ — \X, k) + (X' - X){X' , k'\H\X, k) =0 (A.9) 
'of ' ' 



(A',k'| — |A,k) =0. (A.IO) 



Taking the scalar product of Eq. (|A.7p with |A, k), we obtain 



| = (A..|||A,„)=0. (A.11) 



Since the eigenvalues are time- independent, it is clear that the eigenstates must be time- 
dependent. 

To investigate the connection between the eigenstates of / and the solutions so the 
Schrodinger equation, we first write the equation of motion of |A,k) starting from 
Eq. (|X7ll and using Eq. (|ATT]l : 

(A-/)||A,k) = |^|A,^). (A.12) 

By taking the scalar production with |A', k') and using Eq. ()A.9p to eliminate 

.81 . 

(A',k'||^|A,k) 



we get 

Pi 

\i\ I \ i I 



i{X- X'){X',k'\^\X,k) = {X-X'){X',k'\H\X,k). (A.13) 



^8V 



From this, for A' 7^ A, we infer 



i(A',K'|^|A,K) = {X',k'\H\X,k). (A.14) 
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Eq. (|A.13p does not imply 



i{X',n'\ — \X,n) = {X',k'\H\X,k,). 



If Eq. ()A.14p held for A' = A as well as for A' 7^ A, then we would immediately deduce 
that |A, k) satisfies the Schrodinger equation, that is |A, k) is a special case of 

Note that the phase of |A,k) has not been fixed by our definitions. We are still free 
to multiply | A, k) by an arbitrary time-dependent phase factor. Thus, we can define 
a new set of eigenvectors of I{t) related to our initial set by a time-dependent gauge 
transformation 



where the a\K{t) are arbitrary real functions of time. Because I{t) is assumed not to 
contain time-derivative operators, the |A, orthonormal eigenstates of I{t) just as 

are the |A,k). For A' ^ A, Eq. (|A.13p also holds for matrix elements taken with respect 
to the new eigenstates. Each of the new eigenstates will statisfy the Schrodinger if we 
choose the phases ax^it) such that Eq. (|A.13p holds for A' = A. This requirement is 
equivalent to the following first-order differential equation for the axK{t)- 



Since each of the new set of eigenstates of I{t), | A, k)^, satisfies the Schrodinger equation, 
the general solution is 



where the cak are time-independent coefficients. All of the state vectors with which 
we have dealt so far are time-dependent, while in Eq. ()A.17p we modified the notation 
to indicate the dependence on time explicitly. The Schroodinger state vector is now 
denoted by \t) and the eigenstates of the invariant by \X,K;t). 




(A.15) 




(A.16) 




(A.17) 
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Assume that in the remote past the Hamiltonian H{t) is a constant operator H{—oo) 
having a complete, orthonormal set of time-independent eigenstates \n;i), n being a 
label for all relevant quantum numbers and i standing for "initial state." Similiarly, 
assume that in the remote future, the Hamiltonian is a constant operator H{oo) and it 
possesses time-independent eigenstates |m;/), m labeling the quantum numbers and / 
standing for "final state." The explicit time variation of H{t) for intermediate times is 
arbitrary except for piecewise continuity; in particular, we do not exclude the possibility 
of variations rapid enough to render an analysis in terms of quasistationary states of 
H{t) impossible. 

We want to calculate the transition amplitude T{n — ?• m) connecting an initial state 
\n;i) to a final state |m;/). Thus we consider the case in which the Schrodinger state 
vector I — oo) corresponds to an eigenstate \n;i). The superposition coefficients of 
Eq. (|A.17p for this problem are given by 

CAre = e-^"^'=(-~)(A, k; -oo\n; i) (A.18) 

from which we obtain 

\t) = ^exp(z [axnit) - aAK(-oo)]) |A, k; t)(A, k; -oo|n; z). (A.19) 

The transition amplitude is therefore given by 

T(n — )• m) = (m; /|oo) 

= ^exp(i [ax^ioo) - axn{-oo)]) (m; /|A, k; oo)(A, k; -oo|n;z). (A. 20) 

A, re 

The properties of I{t) apply equally well to any operator that is an invariant corre- 
sponding to a given H{t). In general, for a system of / degrees of freedom, there is 
an infinite family of such invariants, the members of which are functions of a set of / 
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independent invariants. Two such invariants will, in general, have different eigenstates, 
different time derivatives, and different commutators with the Hamiltonian. However, 
one must arrive at the same physical results no matter what invariant we use and, 
therefore, the choice of which particular invariant to use may be made on the basis of 
mathematical convenience. Here we demonstrate that the physical result is independent 
of the choice of invariant, we give a direct proof that a transition amplitude, such as in 
Eq. (|A.20p is indeed independent of our choice of invariant. 

Suppose that we have two complete orthonormal sets of states, \v;t) and \w;t), all 
of which satisfy the time-dependent Schrodinger equation; and suppose that the states 
\v;t) are eigenstates of one set of operators, whose eigenvalues are labeled by v, and 
that the states \w; t) are eigenstates of a different set of operators, whose eigenvalues are 
labeled by w. The transition amplitude T{n — )■ m) can be expressed as 

r(n — )• m) = ^^(m; f\v; oo){v; — oo|n; i) (A. 21) 

V 

or as 

T(n — )• m) = ^^(m; f\w; oo){w; — oo|n; i). (A. 22) 

w 

We want to show directly that these two expressions are the same. The completeness of 
the states \w;t) requires 

\v;t) = ^\w;t){w;t\v;t). (A. 23) 

w 

Operating on this equation with {i{d/dt) — H), and using the facts that all of the states 
satisfy the Schrodinger equation and that the states \w,t) are orthogonal, we obtain 

^^{w;t\v;t) =0. (A.24) 

Thus the quantity {w; t\v; t) is independent of time. We now use the completeness of the 
state \v;t) and \w;t), Eq. ()A.24p . and the orthonormality of the states \w;t) to rewrite 
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Eq.^EMi) as 

T{n^m)= [m] f\'W]oo){w]oo\v-,oo){v-,—oo\'w';—oo){'w'-,—oo\n]i) 

= {m; f\w;oo){w;—oo\v;—oo){v;—oo\w';—oo){w';—oo\n;i) 

= ^^(m; oo)(w; — oo|n; i). (A. 25) 

w 

Thus, Eqs. ()A.2ip and ()A.22p are the same, as asserted. 

Suppose for shnphcity that the eigenstates of / are nondegenerate, so that the eigen- 
value of / is the only quantum number required for describing the system. When this 
is so, as it is in our discussion of the time-dependent harmonic oscillator, then it is 
particularly convenient to choose an invariant having the property that it becomes time- 
independent as t — 7- — oo so that the commutator [/(— oo), oo)] vanishes. Then the 
normalized eigenvectors of H{—oo) and /(— oo) are identical to within constant phase 
factors. Consequently, we may choose the initial state \n;i) simply to be a eigenstate of 
/(— oo), say |A;— oo). Eq. (|A.20p then reduces to 

T{n m) = exp {i [a„(oo) - a„(-oo)]) (m; /|A„; oo), (A.26) 

and the transition probability is given by 

Pnm = \T{n m)|^ 

= Km; /|A„;oo)|2. (A.27) 

As t —7- oo, the invariant operator I{t) in general remains time-dependent and does 
not commute with the Hamiltonian. Therefore, the state |Arra;oo) in Eq. (jA.27p is a 
superposition of eigenstates of H{oo)] this is another expression of the fact that energy 
is not conserved in our system. 
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Prom the structure of Eq. ()A.20p . it is apparent that we may express the transition 
ampUtude as a matrix element of an S matrix by writing 

S = Y^ e*"^^(°°^ |A, k; oo)(A, k; -oole-^^-f— ) , 

T{n^m) = {m;f\S\n;i). (A.28) 
It is easily verified that this operator is unitary: 

S^S = SS^ = 1. (A.29) 

In the special case that the Hamiltonian operators in the remote past and distant future 
are identical, H{—oo) = H{oo), so that the initial and final states are the same set, we 
may define an elastic scattering operator R in the standard fashion: 

5 = 1 + 2TTiR. (A.30) 

The operator R describes the nondiagonal transitions just as S does, but subtracts a 
noninteracting part from the diagonal amplitudes so that {n\R\n) represents a "forward 
reaction amplitude" from the state |n) to the same state. The unitarity of the S matrix 
implies 

y"\{m\R\n)\'^ = -Im{{n\R\n)), (A.31) 

m 

which is a statement of the optical theorem: the total reaction probability is proportional 
to the imaginary part of the forward reaction amplitude. 
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B.l. The Gauss-Codazzi Formalism 

Here we wish to solve Einstein's equations in the presence of stress-energy sources con- 
fined to three-dimensional time-like hypersurfaces for a general metric. Following the 
methods used by Ipser and Sikivie, Ref.|4], we shall use the Gauss-Codazzi formalism. 

The Gauss-Codazzi equations relate the four-dimensional geometry of the overall 
global space-time to their projection onto a three-dimensional hypersurface embedded 
within the original four-dimensional space-time. This is done by investigating the in- 
trinsic and extrinsic curvature of the three-dimensional time-like hypersurface. The 
Gauss-Codazzi formalism allows one to find the equations of motion for a collapsing 
domain wall in a very systematic way. To find the equations of motion, one needs to 
specify the metric (and associated energy-momentum tensor) only. 

In this chapter we wish to develop the Gauss-Codazzi formalism for a general metric 
where the only initial requirement is that the coefficients of the metric depend on position 
and time only. We will then arrive a final equation which depends on the coefficient 
(and derivatives of), as well as its associated energy- momentum tensor, which will allow 
us to find the equations of motion for the collapsing domain wall once the metric is 
completely specified. After we develop the general equations, we will compare our result 
with that found in the literature for two different specified metrics: the Schwarzschild 
and Reissner-Nordstrom metrics, respectively. The Schwarzschild and Reisner-Norstrom 
metric coefficients both depend on position only, hence these are an example of a special 



106 



APPENDIX B. GAUSS CODAZZI 



case of the general method we are working with here. 
B.1.1. The Equations 

Here we fohow the technologies developed in Ref.[l]. Let S denote a three-dimensional 
time-like hypersurface containing stress-energy and let be its unit spacelike normal 
(^(j^" = 1). The three-metric intrinsic to the hypersurface 5 is 

Kb = gag - CaCb (B.l) 

where gab is the four-metric of the space-time. Here hab is known as the projected tensor 
for the hypersurface S, see Ref . |25j . This is due to the fact that, when acting hab on a 
vector v"", it will project it tangent to the hypersurface, hence orthogonal to 

{habV^)e = gabV^e - UbV^e 

= 0. 

Let Vq denote the covariant derivative associated with gab and let 

Da = ha'Vb, (B.2) 

hence Da is the covariant derivative on the induced three-dimensional hypersurface. The 
extrinsic curvature of S, denoted by iTab, is defind by 

T^ab = DaCb = T^ba- (B.3) 

The extrinsic curvature depends on how the hypersurface is embedded in the full four- 
dimensional space-time. The extrinsic curvature is used to differentiate different topolo- 
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gies. For example, intrinsic geometry of a cylinder and a torus can be flat, however, we 
know the exterior geometry of each is different. This different topology is given in the 
extrinsic curvature, which will tell us that we are actually on a torus or a cylinder. 
The contracted forms of the first and second Gauss-Codazzi equations are then given 

by 

^R + TTabTT'^'-Tr'' = -2Gabe^' (B.4) 
habD.TT''^ - DaTT = Gbch\e ■ (B.5) 

Here is the Ricci scalar curvature of the three-geometry hab of S, vr is the trace of 
the extrinsic curvature, and Gj' is the Einstein tensor in four-dimensional space-time. 

Here we will be working with infinitely thin domain walls. The stress-energy tensor 
Tab of four-dimensional space-time then is assumed to have a (5-function singularity on 
S. This in turn implies that the extrinsic curvature has a jump discontinuity across 5, 
since the extrinsic curvature is analogous to the gradient of the Newtonian gravitational 
potential. Therefore we can introduce 

lab = T^+ab - T^-ab (B.6) 

which is the difference between the exterior and interior extrinsic curvatures, and 

Sab = j dlTab, (B.7) 

where / is the proper distance through 5 in the direction of the normal and where the 
subscripts it refer to values just off the surface on the side determined by the direction 
of ib^''. Hence the direction for, say -|-,^'^ will be in the direction of the exterior geometry 
of the domain wall, while — will denote the direction of the interior geometry of 
the domain wall. As we shall discuss below, these geometries will be different for the 
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case of the spherically symmetric domain wall. Using Einstein's and the Gauss-Codazzi 
equations, one can show that (see Ref.[7]) 

Sab = (7ab - hablc") • (B.8) 

We can also introduce the "average" extrinsic curvature 

T^ab = ^ (TT+afe + T^~ab) (B.9) 

which will be important later. 

B.1.2. The Surface Stress- Energy Tensor 

Here we restrict ourselves to sources for which the stress energy tensor is given by, see 
Ref.g] 

5"^* = au^'u^ - r] [h"^ + u'^u''^ (B.IO) 

which is the material sources consisting of a perfect fluid. In Eq. ()B.10p is the four- 
velocity of any observer whose world line lies within S and who sees no energy flux in 
his local frame, and where a is the energy per unit area and r] is the tension measured 
by the observer. For a dust wall it is well known that 77 = 0, while for a domain wall 
rj = a. For a domain wall Eq. (|B.10p reduces to 

5"^ = -ah^K (B.ll) 

We also note that the four-velocity is a time-like unit vector orthogonal to the space- 
like unit normal i.e., 

UaU'' = -l, CaU'' = 0, W = +i- 
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B.1.3. Attractive Energy 

Here we derive equations for an observer who is hovering just above the surface S 
on either side. Let the vector field be extended off S in a smooth fashion. The 
acceleration 

= h\u°-VaU^ -i^U^U^^ab (B.12) 

has a jump discontinuity across S since the extrinsic curvature has such a discontinuity. 
The perpendicular components of the accelerations of observers hovering just off S on 
either side satisfy 



2u''u^^ab 

2l{h''^ + u''u^)^ab-2-S'''>^ab (B.13) 
a a 



and 



-U^U^ab 



AnGn{(T-2r]). (B.14) 



Here we comment on the precense of the second term on the right hand side of Eq. ()B.13p . 
This term takes into account the contributions to the energy-tensor Tab which are present 
in the vacuo on opposite sides of S. For example, if there is only mass present, then Tab 
vanishes off the shell, hence the second term is zero. In the case of charge present, then 
Tab does not vanish, then the contribution to Tab outside can be taken from the Maxwell 
tensor. 
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B.2. Spherical Walls 



In this section we shall obtain the asymptotically flat solutions to Einstein's equations 
for spherically symmetric domain walls with an arbitrary metric. Here we will consider 
two cases. First we will consider the case where the metric coefficients only depend on 
the radial position of the domain wall. Second, we will consider the case where the 
metric coefficients depend on both the radial position of the domain wall and the time. 

B.2.1. Radial dependence only 

For a spherical shell of stress-energy, let the unit normal ^4. point in the outward radial 
direction. It is well known that asymptotic flatness and spherical symmetry requires 
that the interior geometry is flat (Birkhoff's theorem). For the external geometry we 
will choose an arbitrary metric. First we shall consider the case where the coefficients 
only depend on position. Hence, 



A{r)dt^ + B{rfdr^ + r^dfl^ for r > R{t) 



(B.15) 



and 




(B.16) 



where 



dn"^ = de^ + sin^ ed(f>^ . 



(B.17) 



Here the equation of the wall is 



r = R{t). 



(B.18) 
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One finds for tlie components of and (a = t or T, r, 6, cj), in tliat order) 

= (/3^(r)-\ i?,, 0, 0), K_) = (a, i?,, 0, 0), 

= 0, 0), (r-) = «, 0, 0). (B.19) 

Here Rr = dR/dr, wliere r is the propertime of an observer moving with four- velocity 
ti" at the wall, and 



a 



/3 =Atr = ^/Air) + A{r)B{r)R'^. 



(B.20) 
(B.21) 



However, here we should comment that the condition that ^ is of unit normal, this then 
implies the condition that 



B{r) 



A{r)' 



Therefore we can rewrite Eq. ()B.2ip as 



(B.22) 



a 



=Tr = VT+M, 



P = = ^A{r) + ii2. 



(B.23) 
(B.24) 



These expressions and the definitions Eqs. (|B.2p . (|B.3P and ()B.9P imply that 



1 

R' 



(B.25) 



and 



1 
1 

a 



Rr- 
Rtt 



A' 



(B.26) 
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where 

dA{r) 



A' 



dr 



r=Rit) 



(B.27) 



Substituting into Eqs. ()B.13p and (]B.14p then yields the equations of motion 

(a + = - 2^^^^i^ - ^ - (B.28) 

a K la 

aA' 

(a - l3)Rrr =4Tra/3G{a - 2rj) - —. (B.29) 



Taking the ratio of Eqs. (|B.28p and ()B.29p ahows us to ehminate Rrr from the expression, 
so we then find 

{l-A)r] A'a {l-A)S''^nab 

Here we make some general comments on Eqs. (jB.28p and (|B.29p . First, in the absence 
of stress-energy outside the domain wall, R-j-r is always negative provided r/ > 0. Hence 
a spherical domain wall with, say only mass, with 7] > will always collapse to a black 
hole, regardless of its size. Second, in the presence of stress-energy outside the domain 
wall, Rrr is always positive provided that the source term is small compared to the 
other terms. However, if the source term is large compared to the other terms, Rrr can 
become positive at some point. This means that the collapsing object will turn around 
and begin to expand. 

Eq. (]B.30P allows us to find the equations of motion for a specific geometry, provided 
that the coefficient A = A{r), i.e. is only a function of position. In the next section we 
demonstrate the findings in Eq. ()B.30p for two specific cases found in the literature. This 
will allow us to demonstrate ease of the general form of the equations of motion. 
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B.2.2. Radial and Time dependence 



In this section we will write the exterior metric, Eq. (|B.15p . as 



A{r, t)dt^ + B{r, tfdr"^ + r'^dn'^ for r > R{t) 



(B.31) 



where we will maintain that the interior metric is still given by the Minkowski line 
element. We will again take that the equation of the wall is given by Eq. (|B.18p . this 
then gives that the components of and are unchanged in form from Eq. ()B.19p . 
Note however that one does have to make the change from A{r) and B{r) to A{r, t) and 
B{r,t), respectively. As in the case of radial dependence only, the condition that is a 
normalized space-like vector, we again have the condition that 



Therefore we define a and /3 in the same manner as in the case with only radial depen- 
dence, using the suitable substitution. 

We can then find that the acceleration outside and inside the domain wall are given 



B{r,t) 



1 



(B.32) 



Air,t)- 



by 



+ 




A^ + 2R^{A- 13) -?,AI3 



(B.33) 



a 



where A' is given in Eq. (|B.27p and 



A 



dA 
It' 



(B.34) 



114 



APPENDIX B. GAUSS CODAZZI 



Comparing the acceleration outside the domain wah for the radial and time dependent 
metric coefficient, Eq. (|B.33p . to that of the acceleration outside the domain wall for the 
radially dependent metric coefficient, Eq. ()B.26p . we see that the acceleration outside 
the domain in the new case is just the acceleration in the radial case modified by an 
additional term which depends on t-derivatives of the metric coefficient. This is not an 
unexpected result. 

Substituting Eq. (jB.33p into Eqs. (|B.13p and (|B.14p then yields the equations of motion 



(a + l3)Rr 



r]aP{a + P) aA' a/3 . 

2 2 O TTnh 

a R 2 a 

A' + 2R\A-f3)-3Af3 



2A3 L 



(a - p)Rrr =4-^0^0 Nio- - 2l]) 



aA' 



RAa 



A^ + 2R^{A- P) -3AP 



(B.35) 



(B.36) 



Taking the ratio of Eqs. (|B.35p and ()B.36p allows us to eliminate Rrr from the expression, 
so we then find 



=o-(cr - 2r]) + 



(1 - A)7] 



A'a 



+ 



2TT{a + f3)GNR 

(1 - A)S''^^a^ 

ARa 



47r(a + f3)GN 



4TiGNA^{a + P) 



A^ + 2R?{A- f3)-2,Ap 



(B.37) 



Here we make some general comments on Eqs. (|B.35p and (|B.36p . First, we again see 
that the first three terms in Eq. ()B.35p and the first two terms in Eq. ()B.36p are identical 
to the radially dependent metric coefficients only, where the last term comes from the 
time dependence of the metric coefficients. Second, it is not as obvious in this case the 
behavior of the domain wall. In the case of gravitational collapse, < 0, making the 
last term positive. 
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B.3. Examples 

In this section we present some examples using the equation of motion in Eq. ()B.30p . 
First, we will investigate the case of a massive domain wall. We will show that Eq. (|B.30p 
automatically leads to the equation of motion arrived at by Ipser and Sikivie, see Ref.[3]. 
Second, we will investigate the case of a massive-charged domain wall. We will show 
that Eq. (|B.30p automatically leads to the equation of motion arrived at by Lopez, see 
Ref.[6]. 

Here we note that the usual procedure for determining the metric coefficients is to 
consider the asymptotic region of space-time (see for example Ref . |42j ) . Here one writes 
the Ricci tensor, which gives the equations of motion for the the metric coefficients. 
Then using the asymptotic requirements of the space-time, one integrates the equations 
of motion for the metric coefficients and fixes the integration constant. As stated above, 
we will just start with the metric coefficients to find the conserved quantities for the 
collapsing domain wall. 

B.3.1. Massive Domain Wall 

It is well known that asymptotic flatness and spherical symmetry require the exterior 
geometry to be Schwarzschild. Therefore we can write Eq. ijB.lSp . the exterior metric, as 

ids^)+ = -ll —jdt^+il —j dr^ + r^dn^. (B.38) 

Comparing the Schwarzschild metric with Eq. ijB.lSp . one can then identify 

2GM 

A(r) = l . B.39) 

r 

Since the domain wall only contains mass, the stress-energy is only present on the domain 
wall. Hence Tab vanishes outside of the domain wall. Therefore using Eq. ()B.30p we can 
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immediately write 



, 2GM (a - 2r]) 



or rearranging the terms we have 

1 



2GM 



Vi+^+\/i--^ + ^? 



AnaR^ (B.41) 



where in the second hne we use the definition of a and /3, Eqs. (|B.23p and ()B.24p re- 
spectively. This is identical to Eq.(3.8) in Ref. ^j, for the case of the massive domain 
wall. 

B.3.2. Massive- Charged Domain Wall 

Since the domain wall is charged, and spherically symmetric, the geometry outside the 
domain wall is given by the Reissner-Nordstrom solution to Einstein equations. There- 
fore we can write Eq. ()B.15p . the exterior metric, as 

Comparing the Reissner-Nordstrom metric with Eq. ()B.15p . one can then identify 

2GM Q2 

A{r) = 1 — + ^. (B.43) 

In this case the domain wall contains both mass and charge, thus the stress-energy out- 
side of the shell is taken from Maxwell's tensor, since the inside portion of the spherically 
symmetric domain wall will not feel the influence of the charge. The only nonvanishing 
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components outside the domain wall are 



rp 



= Ji 



(B.44) 



By taking the difference of Eq. (|B.4p on opposite sides of S, one finds 



^ oab- _ 

''''' ~ AnaR^- 



(B.45) 



Therefore using Eq. (|B.30p we can write 



{a - 2ri) + 



47r(a + 



27r(a + /3)i?2 



0. 



(B.46) 



Although this is an algebraic equation of second degree in a, only one of the two roots 
holds 

{GM - g) 



a 



27r(a + /3)i?2 ' 

which, using Eqs. (|B.23p and ()B.24p can be put in the form 



(B.47) 



a -(3 = A-naGR. 



(B.48) 



Therefore, solving for the mass yields 



(B.49) 



which is identical to Eq.(61) in Ref.[6], for the case of the massive-charged domain wall. 
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This work was originally completed in Ref . |30| . here we will outline the results. 
In the range t < 0, cj is a constant and the solution to Eq. (|3.5p is 



1 



(C.l) 



In the range of interest, during the time of gravitational collapse, we do not have an 
analytical solution to Eq. ()3.5p . However, we can find certain useful properties of p{t). 
First note that in terms of r] 

(C.2) 



1 - v/Rs ■ 

Then after rescaling, Eq. ()3.5p can be written as 



dr]' 



.12 



f 



(C.3) 



where rj' = rj/Rs, f = ^/oJoP■ The boundary conditions are then 



/(0) = i. ffl = o. 



dr]' 



(C.4) 



The last term in Eq. ()C.3p becomes singular as / — )■ 0. We can then consider a more 
well behaved function for 1//^. For example 



drj'"^ 



-{ujqRs 



1 - ?/ 



(C.5) 
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with boundary conditions 



S(0) = 1, « = 0. (CO) 



Eq. fjC.Sp implies that g{i]') is a monotonicahy decreasing function as long as g{r]') > 0. 
Furthermore, it is decreasing faster than the solution for / as long as / < 1, since the 
1/ in Eq. ()C.3p is a larger "repulsive" force than the g term in Eq. ()C.5p . So 

giv') < firi') (C.7) 

for all r]' such that g{r]') > 0. 

Eq. (]C.5P with initial conditions Eq. ()C.6P can be solved in terms of degenerate hyper- 
geometric functions. The important part for us is that g is positive for all rj' and, in 
particular, ^(1) > for all the values of loqRs that we have checked. Therefore /(ry') is 
positive, at least for a wide range of ojqRs. 

We can find some more properties of p{t). Let /i = /(I) ^ 0. Then the equation for 
/ can be expanded near rj' = 1. 



This shows that 



df ... „ ,2. 



/l 1 
1 - f]' fl 



(C. 



{uoRsYfi ln(l - 7?') ^ -oo (C.9) 
as r/' — )• 1. 

Hence p{r] = Rg) is strictly positive and finite while pr^{r] = Rg) = — oo for finite and 
non-zero dq. Since / = ^/uJop and / — )■ 1 for — ^ 0, we also see that p — )■ oo and — )• 
as ojQ — )■ 0. 

In the range tf < t, u is a constant. However, the solution for p is not constant, unlike 



in the range t <0, since the constant solution l/-y/cJ(tj) does not necessarily match up 
with p{tf—) to ensure a continuous solution. Yet it is easy to check that in this range 
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N = and so there is no change in the occupation numbers. So we need only find 
N{tf — ,uj) to determine N(t — > oo,a)). 
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To get an understanding of the number of particles created in the region near the horizon 
we need to investigate the behavior of the function p{t) near the Schwarzschild radius. 
Near the horizon we can then write the velocity term as 

\Rr\ ^ const = A. (D.l) 
In this limit the position of the shell is then, from Eq. ()5.2ip 

R{t) ^Rq-At (D.2) 
where, as stated in Chapter O Rq is the initial position of the shell, we can write 



\Rt 

Therefore the frequency becomes 



. a (D.3) 



(D.4) 



Therefore the auxiliary equation becomes 
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or using Eq. (j7.37p we can write this as, 



1 dp 2 ^ 

C2 (iT2 ^'^OC((i?o-i?.) 



(D.5) 



After rescaling we can write this as 



where r' = At/{Ro — Rs), and / 
are then 

/(O) 



A^ufR'J'C^/^ 



/ 



1 - t' 



1 

7^ 



(D.6) 



ojq{Rs/C{Rq — Rs)Y/^p. The boundary conditions 

dm 



0. 



(D.7) 



C{Ro-Rs)J ' dr' 

The last term with the 1//^ becomes singular as / — )• 0. Let us consider another equation 
with this term replaced by another more well behaved function. For example consider, 



dT'^~ (Ro-RsV^^ 



f 



1-t' 



(D.8) 



where the boundary conditions in Eq. ()D.7p become 



5(0) 



C(i?o — Rs 



1/4 



dgjO) 
dr' 



0. 



(D.9) 



Eq. ljD.SP implies that g{T') is a monotonically decreasing function as long as ^(t') > 0. 
It is decreasing faster than the solution for / as long as / < 1, since the 1//^ term in 
Eq. (|D.6p . Therefore we have 

fir') > 9{r') (D.IO) 

for all r' such that g{T') > 0. 

The solution for g is positive for all r' and, in particular, g{l) > for all the values 
A^Uq^^C^/^rI^'^ /{Ro- Rs)^'^/^ that we have checked. Therefore /(r') is positive, at least 
for a wide range. 
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Let /i = /(I) / 0. Then the Eq. (|D.6p can be expanded near r' = 1, 



Integrating Eq. ()D.lip we can then write 



/i 1 
1-t' ft 



(D.ll) 



/iln(l -r') ^ -oo (D.12) 



as t' —7- 1. Hence p{t = {Rq — Rs)/A) is strictly positive and finite while Pt{t = 
{Rq — Rs)/A) = —CO for finite and non-zero wq- 

We are calculating the occupation number as a function of frequency uj at some 
fixed time. From Eq. ()D.4p we see that, in order to keep oj fixed in time, ~^ as 

— 7- 0. Thus, ijj varies with ojq and not with time. Since / = {Rs/C'{Rq — 
and / = {Rs/C{Ro — i?s))^/^ for ujq — 0, we see that /? — oo and pr = as — ^ 0. 
This implies taht the occupation number in Eq. (|3.2ip diverges as r —)• Tc since B ^ 
as r — )• Tr- 
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